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Abstract 

We consider the focusing (attractive) nonlinear Schrodinger (NLS) equation with an external, 
symmetric potential which vanishes at infinity and supports a linear bound state. We prove that 
the symmetric, nonlinear ground states must undergo a symmetry breaking bifurcation if the 
potential has a non-degenerate local maxima at zero. Under a generic assumption we show that 
the bifurcation is either subcritical or supercritical pitchfork. In the particular case of double- 
well potentials with large separation, the power of nonlinearity determines the subcritical or 
supercritical character of the bifurcation. The results are obtained from a careful analysis of 
the spectral properties of the ground states at both small and large values for the corresponding 
eigenvalue parameter. 

We employ a novel technique combining concentration-compactness and spectral properties 
of linearized Schrodinger type operators to show that the symmetric ground states can either 
be uniquely continued for the entire interval of the eigenvalue parameter or they undergo a 
symmetry-breaking pitchfork bifurcation due to the second eigenvalue of the linearized operator 
crossing zero. In addition we prove the appropriate scaling for the L q , 2 < q < oo and H 1 norms 
of any stationary states in the limit of large values of the eigenvalue parameter. The scaling 
and our novel technique imply that all ground states at large eigenvalues must be localized near 
a critical point of the potential and bifurcate from the soliton of the focusing NLS equation 
without potential localized at the same point. 

The theoretical results are illustrated numerically for a double-well potential obtained after 
the splitting of a single-well potential. We compare the cases before and after the splitting, 
and numerically investigate bifurcation and stability properties of the ground states which are 
beyond the reach of our theoretical tools. 



1 Introduction 

Over the past few years, there has been a remarkable growth of interest in the study of nonlinear 
Schrodinger (NLS) equations with external potentials. This has been fueled, to a considerable 
extent, by the theoretical and experimental investigation of Bose-Einstein condensates (BECs) 
[261128] . Localized waveforms emerge within the atom-trapping potentials in such ultracold systems 
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[3]. Another major area of applications for NLS equations is nonlinear optics, in particular photonic 
crystals and optical waveguides [2TJ E] • 

One generic type of the external potential for the NLS equation that has drawn considerable 
attention is the symmetric double-well potential. This is due to its relative simplicity which often 
makes it amenable to analytical considerations, but also due to the wealth of phenomenology that 
even such a relatively simple system can offer. Such potentials in the atomic physics setting of BECs 
have already been experimental realized pQ through the combination of routinely available parabolic 
and periodic (optical lattice) potentials. Among the interesting phenomena studied therein were 
Josephson oscillations and tunneling for a small number of atoms, or macroscopic quantum self- 
trapping and an asymmetric partition of the atoms between the wells for sufficiently large numbers 
of atoms. Double well potentials were also examined in the context of nonlinear optics, e.g. in twin- 
core self-guided laser beams in Kerr media [3], optically induced dual-core waveguiding structures 
in a photorefractive crystal |15j , and trapped light beams in a structured annular core of an optical 
fiber [25]. 

In the present work, we address the NLS equation with a symmetric potential as the prototypical 
mathematical model associated with the above experimental settings. For simplicity we will focus 
on the case of one space dimension. We write the equation in the normalized form 

iut = — u xx + V(x)u + a\u\ 2p u, (1.1) 

where u(x, t) : R x R — > C is the wave function, p > is the nonlinearity power, a S R determines 
the defocusing (repulsive), respectively focusing (attractive), character of the nonlinearity when 
a > 0, respectively a < 0, and V{x) : R H > R is an external real-valued, symmetric (even in x) 
potential satisfying: 

(HI) V(x) e L°°(R), 

(H2) lim^i^ V(x) = 0, 

(H3) V(-x) = V(x) for all x € R. 

Hypothesis (HI) implies that —d 2 + V(x) is a self-adjoint operator on L 2 (R) with domain H 2 {M). 
We will make the following spectral assumption: 

(H4) Lo = —d 2 + V(x) has the lowest eigenvalue — Eq < 0. 

It is well known from the Sturm-Liouville theory that all eigenvalues of Lq are simple, the corre- 
sponding eigenfunctions can be chosen to be real valued and the one corresponding to the k + 1-th 
eigenvalue has exactly k zeroes, and, because of the symmetry (H3), is symmetric (even in x) if k is 
even and anti-symmetric if k is odd. We can choose a normalized eigenfunction, ipo, corresponding 
to the eigenvalue — Eq, which will satisfy: 

- Vo(x) + V(x)Mx) + EoMx) = 0, Mx) > 0, M-x) = ifo{x), x G R, ||^o||^ = 1. (1.2) 

We are interested in understanding properties of stationary, symmetric and asymmetric states of 
(jl.ip . i.e. solutions of the form u(t,x) = e lEt (p(x), where (f) satisfies the stationary NLS equation 

- <j)"{x) + V{x)(j){x) + <j\(j){x)\ 2p 4){x) + E<j){x) = 0, xeR, (1.3) 
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and E G R is an arbitrary parameter. We recall the following basic facts about solutions of the 
stationary NLS equation in one dimension. 

(i) Via standard regularity theory, if V(x) G L°°(R), then any weak solution <p(x) G iT (R) of 
the stationary equation (|1.3p belongs to # 2 (M). 

(ii) All solutions of the stationary equation (jl.3p in i7 2 (R) C*(R) are real- valued up to multi- 
plication by e ie , 8 G R. 

(iii) If E > 0, all solutions of (II. 3p in H 2 (M.) decay exponentially fast to zero as \x\ — > oo. 

Numerically we will focus on a one parameter double well potential V = V s constructed from 
splitting the single- well potential Vq(x) = — sech 2 (x) : 

V s (x) = V (x + s) + V (-x + s), s>0. (1.4) 

The general theory of bifurcations from a simple eigenvalue of the linearized operator [23] implies 
that solutions with small H 2 norm of the stationary equation (jl.3p exist for E near Eq. The 
symmetry hypothesis (H3) implies that these solutions are symmetric (even in x). Variants of the 
local bifurcation analysis near E = Eq, including the fact that E > Eq if a < 0, and E < Eq if 
a > 0, have already appeared in [3T , 27J, as well as in many recent publications. We review this 
analysis in Section 2 to give readers a complete picture. 

Orbital stability [32] of the stationary state e iEt <j)(x), (j){x) G R, is closely related to the lin- 
earization of the time-dependent NLS equation (jl.ip at the stationary state, which, in the direction 
e lEt [ui(x,t) + ivQ(x,t)], is given by: 
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where L± are self adjoint linear Schrodinger operators with domains H 2 {W) C L 2 (R) : 
/ L + (<p,E) = -d 2 + E + V(x) + a(2p + l)\<p\ 2 P(x), 

\ L_(cP,E) = -d 2 + E + V(x) + a\cP\ 2 P(x). U ' 0j 

Sufficient conditions for orbital stability and orbital instability, which we will use throughout this 
paper, were obtained in [32], [9] [8]. 

Definition 1 . If (<f>,E) solves (11. 3p and zero is the lowest eigenvalue of L_ ((f), E) we call <f> a ground 
state of (fOp . 

Remark 1. Note that, for any solution (<p,E) of (jl.3p . zero is an eigenvalue of L^((p,E) with 
eigenfunction 4>. Via standard theory of second order elliptic operators the above definition is equiv- 
alent to the one requiring a ground state to be strictly positive or strictly negative. 

In particular it is known that solutions (<p } E) of the stationary NLS equation (jl.3p with \E — Eq\ 
and II^Hh 1 small are orbitally stable ground states, see Section 2. We remark as a side note that, 
for critical and supercritical nonlinearities, p > 2, and more restrictive hypotheses on the potential 
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V(x), one can show asymptotic stability of these solutions in the space of one dimension jUJ 123] . 
For subcritical nonlinearities asymptotic stability is proven only in dimensions higher than one, see 

nana us). 

Kirr et al. |16j showed that the symmetric ground states undergo a symmetry-breaking bifur- 
cation at E = E* > Eq, in the focusing case a < 0, with p = 1, or other cubic like nonlinearities, 
provided the first two eigenvalues of Lq = —d^ + V(x) are sufficiently close to each other. In partic- 
ular, the result is applicable to double-well potentials such as (jl.4p for sufficiently large separation 
parameter s between the two wells. Furthermore the authors show that the symmetric states be- 
come unstable for E > E*, and, a new pair of orbitally stable, asymmetric ground states exist for 
E > E*. The proofs rely on a Lyapunov-Schmidt type projection onto the two eigenvectors corre- 
sponding to the lowest eigenvalues of Lq, which exists for small ||0||#i, combined with a normal 
form analysis of the reduced system valid to all orders. Marzuola & Weinstein [22\ used a time 
dependent normal form valid for finite time to extract interesting properties of the dynamics of 
solutions of the NLS equation (jl.ip near the bifurcation point E = E*. We also mention that |30| 
uses a similar finite time, normal form technique to study the solutions and predict bifurcations 
of the first excited (anti-symmetric) state for defocusing NLS equation (a > 0) with a symmetric 
double well potential which is essentially brought in the large separation regime by passing to the 
semi-classical limit and assuming a \ at a specific rate. 

To our knowledge there are very few results for bifurcations of NLS stationary states in non- 
perturbative regimes. Rose & Weinstein |29] use variational methods to show that the stationary 
NLS equation (|1.3p with a < 0, and potential satisfying (HI), (H2) and (H4), has at least one 
solution for any E > Eq. Jeanjean & Stuart [TT] prove that for a < 0, the symmetric states 
bifurcating from the lowest eigenvalue of Lq can be uniquely continued for all E > Eq, hence there 
are no bifurcations along this branch, provided V(x) is monotonically increasing for x > 0, and C 1 
in addition to satisfying (H1)-(H4). In particular, the result applies to the potential (jl.4p if s < s*, 
where ^ 

s* = -arccosh(2) = arccosh(v / 3/v / 2) ~ 0.6585, (1.6) 

because for x > and < s < s* : 

V' s {x) = Vq(x -s) + Vq(x + s) = smh(2x) + cQsh ^ cosh ( 2s ) - cosh 2 (2s)] > 0. 

cosh (x — s) cosh (x + s) 

Results on continuation of branches of stationary states in the defocusing case a > 0, but without 
reference to existence or non-existence of bifurcation points can be found in |12|. 113] . In [2] the 
authors rely on variational techniques to deduce symmetry-breaking of the ground states in Hartree 
equations. Their method can be adapted to our problem and implies the emergence of asymmetric, 
ground state branches in the focusing case a < provided the nonlinearity is subcritical, p < 2, 
and V(x) is continuous, bounded, and has at least two separated minima. In particular, assuming 
p < 2, asymmetric ground state branches will appear for the potential (jl.4p as soon as it becomes 
a double well, i.e., for s > s*, but the method cannot tell whether the asymmetric branches are 
connected to the symmetric branch of ground states bifurcating from the lowest eigenvalue of Lq. 
Jackson & Weinstein [10] use a topological shooting method for the case p = 1 and Dirac type 
double- well potential, i.e. Vq(x) = 5{x) in (jl.4p . to show that the asymmetric branches emerge 
from the symmetric ones via a pitchfork bifurcation and they all coexists past a certain value of E. 
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Our main result extends the ones in [IB] to non-perturbative regimes and the ones in [2] to critical 
and supercritical nonlinear ities, p > 2, while proving that the asymmetric ground states emerge 
from the symmetric ones via a pitchfork bifurcation. The main theorem is formulated as follows. 

Theorem 1. Consider the stationary NLS equation fjl .3j) in the focusing case a < 0, with V(x) 
satisfying (H1)-(H4). Then the C 1 curve E i— Y <f> = ipE £ H 2 , E > E$, iPe{x) > of symmetric, 
real valued solutions bifurcating from the zero solution at E = Eq, undergoes another bifurcation at 
a finite E = E* > Eq, ipE* £ H 2 , provided V(x) has a non-degenerate maxima at x = 0, and 

(H5) xV'(x) € L°°(R). 

The bifurcation is due to the second eigenvalue X(E) of L + (ipE,E) crossing zero at E = E*. 
Moreover, if p > 1/2, and the following non- degeneracy condition holds: 

then the bifurcation is of pitchfork type: the set of real valued solutions (</>, E) S H 2 x R in a 
neighborhood of (ipE. t ,Ez) consists of exactly two orthogonal C 1 curves: the symmetric branch 
E i—T- <p = ipE which continues past E = E*, but becomes orbitally unstable, and an asymmetric 
branch (<f)(a), E{a)) : 

E(a) = E* + ^a 2 + o(a 2 ), 0(a) = ipE* + «0* + 0(a 2 ), a £ R, \a\ small 

where <f>* is the eigenfunction corresponding to the second eigenvalue of L + (ipE t , E*), and Q can be 
calculated from tpE, and see (|3.35p . The asymmetric solutions are orbitally stable if Q > and 
||0(a) ||^2 is increasing as E{a) increases, but they are orbitally unstable if [|0(a)||j^a is decreasing 
with E(a) or if Q < 0. 

In particular, the result applies to the potential (|1.4p for s > s* = arccosh(y / 3/v / 2) because 

y s '(0) = Vq(s) - Vq{s) = 0, V"(0) = 2F "(s) = 12sech 4 (s) - 8sech 2 (s) < 0, s > s*, 

and implies that a pitchfork bifurcation occurs along the branch of symmetric states. Recall that 
for s < this branch can be uniquely continued for all E > Eq due to the result in [11] , see (jl.6p . 
Moreover, in the large separation limit s — > oo, the branch of asymmetric states is orbitally stable 
near the pitchfork bifurcation if p < p* , where 

P * = 2 ( 3 + ^) ~ 3 ' 3028 ' ( L7 ) 

and orbitally unstable for p > , see Corollary [23 The threshold power of the nonlinearity was 
predicted in |30j but we justify this result with rigorous analysis. 

We emphasize that hypotheses (HI) and (H5) can be relaxed to V(x), xV'(x) G L <? (R) for some 
q > 1, at the expense of slightly complicating the proofs in this paper. Moreover, our results 
extend to more than one dimension x S 1", n > 2, and other symmetries in R n , provided that 
the symmetries still prevent the solutions to concentrate at |x| = oo, see Remark [71 Note that 
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for classifying the bifurcation, we will have to assume that the second eigenvalue of L + (ipE, E) is 
simple. To completely remove any symmetry assumptions, or the spectral assumption (H4), or the 
simplicity of the second eigenvalue of L + , is a much more difficult task, see [20J for partial results. 

The proof of the main result relies on Theorems [21 El El and [6] which, viewed individually, are 
important themselves. Properly generalized they could completely describe the set of all solutions 
of the stationary NLS equation (|1.3|) for any dimension, arbitrary potentials and more general 
nonlinearities in terms of the critical points of the potential. In Section [2] we prove the following 
dichotomy: the branch of stationary solutions (ipE,E) bifurcating from the lowest eigenvalue — Eq 
of —d 2 + V(x) can either (a) be uniquely continued for all E > Eq or (b) there exists a finite 
> Eq such that zero is an accumulation point for the discrete spectrum of the linearized operator 
L + (ipE, E), as E f - E*. The result essentially eliminates the possibility that ||^.e||#2 diverges to 
infinity as E approaches E* with ||LT 1 (^e, -E)||l 2 i-s>l 2 remaining uniformly bounded, and relies on 
the differential estimates for the mass (charge) and energy of the NLS equation (jl.ip . 

In Section [3] we use a novel technique combining concentration-compactness, see for example 
[5], and the spectral properties of the linearized operator L+, to show that, in case (b), the states 
ipE must converge in H 2 to a nonzero state tpE, for E = E*. By continuity we deduce that the 
linearized operator L + (xpE t , E*) has zero as a simple eigenvalue, then we use a Lyapunov-Schmidt 
decomposition and the Morse Lemma, see for example |24| . to show that a pitchfork bifurcation 
occurs at (ipE,, E*). The symmetry hypothesis (H3) implies that iPe(x) are even in x which is 
essential in showing that the limit lim^ /•e^iPe exists. It turns out that without assuming (H3) ipE 
may drift to infinity as E approaches E*, i.e. there exists ue G M, lim^ /-e,Ve = ±oo, such that 
lim.E y^E, iPe{' — Ue) = ipE, where ipE, is now a solution of the stationary NLS equation (II. 3j) with 
V(x) = and E = E*, see [20] for a more detailed discussion of this phenomenon. 

In Section [H we obtain new, rigorous results on the behavior of all stationary states in the 
focusing case a < 0, for large E. In Theorem[5]we combine Pohozaev type identities with differential 
estimates for mass, energy and the L 2p+2 norm of the stationary solutions to prove how the relevant 
norms of the solutions scale with E as E — > oo. In particular, we obtain the behavior of the L 2 
norm of the symmetric states, which was numerically and heuristically predicted in [29]. Moreover, 
by combining these estimates with our novel concentration-compactness/spectral technique we 
show that, modulo a re-scaling, the symmetric branches of solutions, along which L + has only one 
negative eigenvalue, converge to a non-trivial solution of the constant-coefficient stationary NLS 
equation, 

-(p"{x)+o\(P(x)\ 2p (p{x)+(p{x) = 0, cf)£H 2 (R). (1.8) 

Since the set of solutions of the latter in dimension one is well known, we adapt and extend the 
bifurcation analysis of Floer & Weinstein [7J to our problem and obtain detailed information on all 
stationary solutions, which, modulo a re-scaling, bifurcate from a non-trivial solution of equation 
(jl.8p . We show that such solutions can only be localized near a critical point of the potential V(x), 
and they are always orbitally unstable for supercritical nonlinearities, p > 2. If the critical point is 
a non-degenerate minimum, respectively maximum, then there is exactly one branch of solutions 
localized at that point, and these solutions are orbitally stable if and only if we are in the critical and 
subcritical regimes, p < 2. We note that compared to the semi-classical analysis in [7], completed 
with orbital stability analysis in [91 Example C], we are forced to make a precise analysis up to 
order four, instead of two, in the relevant small parameter. In addition, we prove non-existence 
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of solutions localized near regular points of the potential V(x), uniqueness of solutions localized 
near non-degenerate minima and maxima, and we recover the stability for the critical nonlinearity 
p = 2. All these results can now be extended without modifications to the problem studied in [TJ. 

In Section [5j we illustrate the main theoretical results numerically for the potential (|1.4p , sub- 
critical nonlinearity p = 1 and supercritical nonlinearities p = 3, p = 5 in the focusing case a < 0. 
We note that p = 3 < p* and p = 5 > p*, where p* is defined by (|1.7|) . We will show that both 
subcritical and supercritical pitchfork bifurcations occur depending on the value of p. For this po- 
tential we will also show numerically that, except for the bifurcation predicted by our main result, 
there are no other bifurcations along any of the ground state branches, a result beyond the grasp 
of our current theoretical techniques. 

In what follows, we shall use notations 0(e) and o(e) as e — > in the sense 

A = Oie) ^ lim e -1 A(e) =A oa eR and A = o(e) ^ lim e^Aie) = 0. 

We will denote e-independent constants by C, which may change from one line to another line. We 
will also use the standard Hilbert space L 2 (R) of the real valued square integrable functions on a 
real line and the Sobolev space H 2 (R) C L 2 (IR) of the real valued functions on R which are square 
integrable together with their first and second order weak derivatives. 

Acknowledgments. PGK is partially supported by NSF-DMS-0349023 (CAREER), NSF- 
DMS-0806762 and the Alexander-von-Humboldt Foundation. EWK is partially supported by NSF- 
DMS-0707800. DEP is partially supported by the NSERC. The authors are grateful to V. Natarajan 
and M.I. Weinstein for fruitful discussions, as well as to C. Wang for assistance with some of the 
numerical computations. 

2 Local bifurcations of symmetric ground states 

In this section we trace the manifold of symmetric ground states of the stationary problem (jl.3p 
from its local bifurcation from the linear eigenmode of Lq = —d 2 + V(x) near E = Eq up to its 
next bifurcation. We will show that the symmetric state exists in an interval to the right of Eq 
if a < and in an interval to the left of Eq if a > 0. In the case of a < 0, we will further find 
necessary and sufficient conditions for the symmetric state to be extended for all values of E > Eq 
or suffer a symmetry-breaking bifurcation. 

Let us rewrite the stationary equation (|1 .3j) for real-valued solutions <p{x) as the root-finding 
equation for the functional F(<p, E) : H 2 (R) xl4 L 2 (R) given by 

F(4>, E) = {-d 2 x + V{x) + E)<f> + a^l 2 ^. (2.1) 

We recall the following result describing the existence of symmetric ground states near i?o- 

Proposition 1. Let —Eq < be the smallest eigenvalue of Lq = —d 2 + V(x). There exist e > 
and 5 > such that for each E on the interval I e , I e = (Eq — e, Eq) for a > 0, I e = (Eq,Eq + e) for 
a < 0, the stationary equation <\1.3h has exactly two nonzero, real valued solutions ±iPe( x ) G H 2 (R), 
satisfying \\iPe\\h 3 < ^- Moreover 

U E \\m <C\E-E \k. 
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for some C > 0, the map E i-> ^ is C 1 from X e to H 2 and ?Pe(x) = iPe(-x) f° r each i£R and 

E el £ . 

Proof. We sketch the main steps. For any p > 0, the functional F : H 2 (K) xlh) L 2 (M) is C 1 , 
i.e. it is continuous with continuous Frechet derivative 

D 4> F{4>, E) = -d 2 x + V(x) + E+(2p + l)a\<j)\ 2p . 

We note that D^F(0, E) = Lq + E. Let tpQ be the L 2 -normalized eigenfunction of D^Ffo, E) = 
Lq + Eq corresponding to its zero eigenvalue. Then D^Ffi, Eq) is a Fredholm operator of index 
zero with Ker(Lo + Eq) = spanj-^o} and Ran(Lo + Eq) = [Ker(Lo + Eq)^. Let 

P\\4> = (ifj ,ct>} L 2i> , p±4> = 4>- (iPo,4>)l^o 

be the two orthogonal projections associated to this Lyapunov- Schmidt decomposition. Then 
F(cj), E) = with eft = aipQ + ip, where a = (ipQ, </>) and (p = P±4>, is equivalent to two equations 

P ± (L + E)P ± v + orP±\a4>o + p\ 2p {a-ip + (p) = (2.2) 
(E - E )a + o-(4>o,\a4>o + <p\ 2p (a4> + v)) = (2.3) 

Since P±(Lq + E)P± has a bounded inverse for any E near Eq, the Implicit Function Theorem 
states that there exists a unique C 1 map IR 2 B (a, E) \- > ip = Lp*(a,E) £ H 2 for sufficiently small 
\a\ and \E — Eq\ such that if solves equation (12. 2j) and 

3C>0: \\Ma,E)\\ H 2 <C\a\ 2p+1 . 

Hence (12. 3D becomes a scalar equation in variables (a, E) £ IR 2 given by 

a{E-E ) + a(^ ,\a^o + <f\ 2p (ai> + <p)) L 2 = 0. (2.4) 

Dividing (|2.4p by a and invoking again the Implicit Function Theorem for functions, we obtain the 
existence of a unique continuous map IR B a \— > E = E*(a) £ ]R for sufficiently small \a\ such that 
E solves (H3]) and 

3C> : |£*(a) - S + ^||^o||^+ 2 |o| 2p | < C|a| 4 P. 

Therefore, 2? > i?o if °~ < and E < Eq if <r > 0. Moreover, the map a \— > E*(a) is C 1 and invertible 
for a > 0, rendering the map E 1 t-t ipE '■= aipo + E*{a)) to be C 1 . The negative branch — i/je 
is obtained from the negative values of a. Moreover, since for any solution i/je(x) of (II. 3D we have 
that iPe{—x) is also a solution, uniqueness and continuity in a imply that tPe(-x) = iPe(x). □ 

Remark 2. A result similar to Proposition [7] can be obtained for the anti- symmetric state of the 
stationary equation ( ti.3j) . which bifurcates from the second eigenvalue —E\ of Lq = —d 2 + V(x), 
if the second eigenvalue of Lq exists. Moreover, (0,E) £ H 2 (M.) x R is the only solution of the 
stationary NLS equation (|1.3j) in a small neighborhood of (0, E*) £ H 2 (W) x IR /or any E* > 0, i/ 
—E* is not an eigenvalue of Lq. 
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Let us introduce operators L + and L_ along the branch of symmetric states (ipE,E) according 
to definition (j 1 . 5 1) with <j) = ipE- Since they depend C 1 on E G X e and continuously at Eq, their 
isolated eigenvalues depend C 1 on E G X e and continuously at Eg. 

Since L-^e = and L_ = Lq + Eq for E = Eq, is the lowest eigenvalue of L_ for all E £ I t . 
On the other hand, we have 

L + — L_ = a2pTp^. 

Hence, L + < L_ if a < 0, while L + > L_ if cr > 0, and, via eigenvalue comparison principle, the 
lowest eigenvalue of L + is strictly negative if cr < and strictly positive if a > 0. Consequently, 
is not in the discrete spectrum of L + nor in the essential spectrum for E G X € . The latter follows 
from V(a?) G L°°(R), limi^^ V(s) = 0, and \^ e \ 2 p G L 2 (M), since G fl" 2 (M) L«(R) for any 
2 < c/ < oo. Together they imply that L + — E is a relatively compact perturbation of — <9 2 , hence, 
via Weyl's theorem, the essential spectrum of L + is the [E, oo) interval. 

The following result shows that we can continue the branch of symmetric states (ipE,E) as long 
as is not in the spectrum of L + . 

Lemma 1. Let ip El (x) G H 2 (M) be a real valued solution of the stationary equation (jl.3p for 
E = E\ and assume is not in the spectrum of L + (tpE 1 , E±). Then there exist e > and 8 > such 
that for each E G (E% — e, E\ + e) the stationary equation (II. 3p has a unique, real valued, nonzero 
solution tPe( x ) G iX 2 (1R) satisfying 

Moreover the map E i-> ^ is C 1 /roro (Ei — e, £i + e) to H 2 . 

Proof. The result follows from the Implicit Function Theorem for F((p,E) = at (ip El ,Ei). □ 

Combining the Proposition and Lemma, we get the following maximal result for the branch of 
symmetric modes bifurcating from the point Eq. For definiteness, we state and prove this result 
only for the focusing case a < 0. 

Theorem 2. If a < 0, then the branch of solutions (ipE,E) of (11.31) which bifurcates from the 
lowest eigenvalue —Eq of Lq can be uniquely continued to a maximal interval (Eq,E*) such that 
either: 

(a) E* = oo; 
or 

(b) E* < oo and there exists a sequence {E n } n ^ C (Eq,E*) such that lim n _ >0O £' n = E* and 
L+(jpE„i E n ) has an eigenvalue X n satisfying lim n _ !>00 

Proof. Define: 

E* = sup{E : E > Eq, E h- > ip E is a C 1 extension on (Eq,E) of the map in Proposition [H 
for which is not in the spectrum of L + } 
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Proposition Q] and the discussion following it guarantees that the set above is not empty. Assume 
neither (a) nor (b) hold for E* > Eq defined above. Then we can fix any E\, Eq < E\ < E* and 
find that, for E G [E\, E*), the spectrum of L + (vp E , E) (which is real valued since L + is self-adjoint) 
has no points in the interval [—d,d] for some < d < E\. Indeed, as discussed after Proposition 
CD the essential spectrum of L + at E is [E,oo), and if no d > exists, there must be a sequence 
of eigenvalues X n for L + at E n G [E\,E*) such that lim n _ 5 . 00 A n = 0. But [E\,E*\ is compact 
since (a) does not hold, hence there exists a subsequence E nk of E n converging to E 2 G [Ei,E*\. 
However, E 2 / E* because we assumed that (b) does not hold, so E 2 G [E%,E^) and by continuous 
dependence of the eigenvalues of L + on E G [E\ , E*) we get that is an eigenvalue of L + at 
E 2 < E* which contradicts the choice of E*. 

Consequently L^ 1 : L 2 {U.) h-> L 2 (M) is bounded with uniform bound K = 1/d on [Ei,E*). 
Moreover, by differentiating (jl.3j) with respect to E we have: 

L + d E ^E = -^E d E ip E = -L+^je, Ee(E ,E*), (2.5) 

hence 

WOe^eWl^ < KUeWl^, Ee[E 1: E*). (2.6) 
and, by Cauchy-Schwarz inequality: 

^Ue\\ 2 L 2 = 2{d E il>E,il>E) < 2KU E \\ 2 L2 , E G [Ei,E # ). 

The latter implies 

U E \\h < U El \\W K(E - El) , E G [E U E*), (2.7) 

which combined with E* < oo and bound (I2.6P gives that both Oe^e and ipE have uniformly 
bounded L 2 norms on [Ei,E*). By the Mean Value Theorem there exists ?Pe*(z) G L 2 (R) such 
that 

lim Ue-^eJl^ =0. (2.8) 

We claim that ^^(R) G i/ 1 (IR) is a weak solution of the stationary equation (jl.3jl with E = E*, 
hence it is in .ff 2 (R) since V G L°°(R). Indeed, consider the energy functional 

£{E)= [ \Vi, E {x)\ 2 dx+ [ V( x )\^ E (x)\ 2 dx + ^— [ \i, E {x)\ 2p+2 dx. (2.9) 

JR JR P + 1 JR 

Note that because, ip E is a weak solution of the stationary equation fjl .3 j) for any S G (i^O) ^*) we 
have 

-^ = -1E($eMe)i? (2.10) 

and, via Cauchy-Schwarz inequality: 

d£ 



d£ 



<EU E \\ L ,\\d E ^ 



e\\l 2 - 



Hence, from the uniform bounds (|2.6|) and (|2.7p . the derivative of £(E), and hence £(E), is uniformly 
bounded on [£i , E 1 * ) . On the other hand, from the weak formulation of solutions of (11.31) . we get 

IIW> E ||| 2 + I V{x)\ii E {x)\ 2 dx + a\\^ E \\% + P l 2 + E\\i, E \\ 2 L2 = (2.11) 

JR 
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Subtracting the latter from (p + 1)£ we get that there exists an M > such that 

p\\V^ E \\h +P f V(x)\i> E {x)\ 2 dx <M, for all Ee [Eh. ,£*). (2.12) 



From Holder inequality we obtain: 



V(x)\il> E {x)\ l dx 



< ||^I|l°°IIV ; eIIl 2 - 



Using the inequality in (|2.12p we deduce that ||VV>.e||l 2 has to be uniformly bounded. Consequently, 
there exists M > such that 

H^bIIh 1 < M, for all E G [E^E*). (2.13) 
Because of the embedding of H l {M) into L°°(R) and the interpolation 

< llfll 2 /, 5 llfllr" 2/5 , a>2, 



li ^ \\j\\ L 2 imi L oo' , 



bound (|2.13p together with convergence (|2.8p imply that as E /* E* we have: 
ipE — > i>E t , in L 2 CR) } ii- ro/nN 



Now, by passing to the limit in the weak formulation of the stationary equation (jl.3p . we conclude 
that ipE* ( x ) G H 1 (M) is a weak solution. Moreover, the linearized operator L + depends continuously 
on E on the interval [E\,E*]. By the standard perturbation theory, the discrete spectrum of L + 
depends continuously on E. Since is not in the spectrum of L + for E G \E\, E*) and we assumed 
that (b) does not hold we deduce that is not an eigenvalue of L + at E*. Moreover, since the 
essential spectrum of L + is [E, oo), is not in the spectrum of L + at E = E*. Applying now 
Lemma [1] we can continue the C 1 branch (ip E ,E) past E = E* which contradicts the choice of E*. 

The theorem is now completely proven. □ 

Remark 3. Let a < and L + be computed at the branch points {ipE,E) for E G (Eq,E*), where 
E* is given by Theorem^ Then, L + has exactly one (strictly) negative eigenvalue. This follows 
from the fact that the eigenvalues of L + depend C 1 on E G (Eq,E*), is not in the spectrum of 
L + and for E near Eq, L + has exactly one strictly negative eigenvalue, see the discussion after 
Proposition d 



3 Symmetry-breaking transitions to asymmetric states 

In this section we show that, in the focusing case a < 0, the second alternative in Theorem [2] 
occurs if and only if the second eigenvalue of L + , \(E) crosses the zero value at E = Moreover, 
at E = E*, under the generic assumption of \'(E*) / 0, the branch of symmetric states suffers a 
symmetry-breaking bifurcation of pitchfork type with a new branch of asymmetric states emerging. 
The new branch consists of solutions of (jl.3p that are neither even nor odd in x. Depending on the 
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sign of a quantities Q and R, see Theorem [U which can be numerically computed, the asymmetric 
solutions are either orbitally stable or orbitally unstable with respect to the full dynamics of the 
NLS equation (jl.ip . For double well potentials with large separation, e.g. (|1.4|) with s — > oo, 
the orbital stability of the asymmetric branch is determined by the power of the nonlinearity, see 
Corollary [2J Regarding the branch of symmetric states, it continues past the bifurcation point but 
it is always orbitally unstable. 

Note that sufficient conditions for the symmetry-breaking bifurcation at E* were presented by 
Kirr et al. |16| under the assumption that E* was sufficiently close to Eq and H^H/fi was small 
for all E £ (Eq,E*). Moreover, only the case of cubic nonlinearity (p = 1) was considered. We now 
present a generalization of this result which allows large values of | E* — Eq | , large norms || 'fell if 1 , 
and any power nonlinearities p > 1/2. 

Remark [3] showed that, in the focusing case a < 0, the operator L + on the branch of symmetric 
states has exactly one negative eigenvalue. We first show that this eigenvalue cannot approach zero 
as E approaches E*. 

Lemma 2. Let a < and consider the branch of solutions (fe,-E) of fjl .31) which bifurcates from 
the lowest eigenvalue —Eq of Lq. Let (Eq,E*) be the maximal interval on which this branch can 
be uniquely continued. Lf Xq(E) is the lowest eigenvalue of L + along this branch then there exist 
5, d > such that: 

^o(E) < -d < 0, for all E G [E* - 5, E*) if E* < oo, or for all E > E , if E* = oo. 



Proof. Assume the contrary, that there exists a sequence E n /• E* such that Ao(-E ra ) f~ 0. From 
the min-max principle we have for all integers n : 

X (E n ) = inf I, 1 2 (Tp En ,L + i> En ) 

fcH*, \\t\\ L 2=l WEnWtf 

1 II / II 2 P+ 2 

= 7rr^r&E n ,L^ En + 2p alfej^fej = 2p a ™^ +2 < 0, 

where we used the definitions (|1.5p and L_-fe = for all E £ (Eq, E*). From the above inequalities 
and Ao(-En) /* we conclude that: 

II / l|2p+2 

iim UEjL2 r = o- 



\L 2 



But (ip En ,E n ) solves the stationary NLS equation f)l .3[) . Plugging in and taking the L 2 scalar 
product with .. ^ E ^ 2 we get: 

\\WE n \\ L 2 

^fejggg _ I fe n ( R2 | VM> fen \ „ 
||fej|£ 2 \WE n \\L* H*feJIW 

where the last inequality follows from — Eq being the lowest eigenvalue of Lq = —d 2 + V(x). Passing 
to the limit when n — > oo we get the contradiction Eq > E*. □ 
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Lemma[2] combined with concentration compactness method and the spectral theory of Schrodinger 
type operators enables us to deduce the following important result regarding the behavior of ipE 
for E near E*. 

Theorem 3. Let a < and consider the branch of solutions (ipE, E) of H 1 . 3[> which bifurcates from 
the lowest eigenvalue —E$ of Lq. Let (Eq,E*) be the maximal interval on which this branch can be 
uniquely continued. Denote 

N(E) = |KM£ 2 , Ee(E ,E*). 

If E* < oo, then: 

(i) N(E) is bounded on \Eq,E*), N* := Iiioie^e, N(E) exists, and < iV* < oo; 
(ii) there exists ipE* G -ff 2 (M) such that (tpE*, E m ) solves the stationary NLS equation (II. 3p and: 

lim \\vpE ~ iPe*\\h 2 = 0- 

E->Et 

Proof. For part (i) denote £ E G H 2 (M>), ||^||x 2 = 1 the normalized eigenfunction of L+ corre- 
sponding to its lowest eigenvalue, \§{Ei). We use the orthogonal decomposition: 



i>E = *I>e) He + ^e, where (£° , ^) = 0, E G (E , E* 



Then 



< KOMP < N(E) 
- -Ao(^) --Ao(^)' 

where we used (12, 5p . \q(E) < 0, and the fact that L + > on the orthogonal complement of £ E , 
see Remark [3j The above inequality implies that 

rE dE 

N(E) < N{E l )e iE ^ - A o(B) , for any E < E\ < E < E*. 

Using now the bound < ^, E G [E* — 5,E*) given by Lemma [21 together with the fact that 
N(E) is continuous on [Eq,E* — 5], see Proposition Q] and Theorem [21 we deduce that there exists 
N > such that: 

< N(E) < N, for all E G [E , 



To show that N(E) converges actually to a finite value as E /* E* we go back to the bound ()3.1 
and integrate it from E\ = E* — 5 to any E, E\ < E < E*. We get 



j\^L^)dE = J* l %^y dE + l -N{ El ) - \n{E) (3.2) 



N 1 

< 1[ S + -N(E i: 
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But the integrand on the left hand side is non-negative, hence the uniform bound implies that the 
integral on the left hand side converges as E /~ E*. Since the same holds for the integral on the 
right hand side, we deduce from (|3.2|) that N(E) must converge to a finite limit as E 

Moreover, since both integrals in (|3.2p are now convergent on [E\, E*) we deduce that the deriva- 
tive of N(E) is absolutely convergent on the same interval, see (|3,ip . consequently the derivative 
of the energy functional, see (|2.9|) and (|2.1(jp . is absolutely convergent and the energy functional 
remains uniformly bounded on [Ei, E*). By repeating the argument in the proof of Theorem [2] we 
get uniform bounds in H 1 norm, see (|2.13p . i.e. there exists M > such that: 

UE\\m < M, for all E E [Ei,E*). (3.3) 

Now, if YunE-^E* N{E) = 0, then, equivalently, lims-^ HV'-eIIl 2 = 0- Because of the bound (|3.3p . 
Sobolev imbedding and interpolation in L q spaces, we get lim^^^ HV'eIIl'j =0, q > 2. Hence L + 
depends continuously on E E \E\, E*], and L + at E* becomes Lq + E*. Then the eigenvalues of L+ 
should converge to the eigenvalues of Lq + E* as E /£»,so0 is an eigenvalue of Lq + E* , because 
we are in case (b) of Theorem [2j Equivalently — E* < —Eq is an eigenvalue of Lq, which leads to a 
contradiction. 

Part (i) is now completely proven. For part (ii) we will use the following lemma. 

Lemma 3. Under the assumptions of Theorem^ if there exist a sequence on the branch (^E n ,E n ) 
and a function ip E , E such that 

E n / E, (3.4) 

ipE n ipE, (3.5) 

"4>E n ->■ ipE*, for all 2 < q < oo, (3.6) 

then ipE n 4>E*- Moreover, ipE* is a solution of (II. 3D and: 

lim \\ipE ~ ipE*\\w = 0. 
E^fE, 



Proof. The equation (|1.3p satisfied by the elements on the branch can be rewritten: 

^E n = {-d 2 x + E n )- l [-V{x)i, En {x) - a\^E n (x)\ 2 ^ En ], (3.7) 
and, from (I3.6p . we have 

HE n {x)\ 2 ^ E A\^E t {x)\ 2P ^. 

To show that V(x)ipE n (x) also converges in L 2 (IR) we use the following compactness argument. Fix 
e > and choose R > sufficiently large such that 

11^0*0 IU~ ({!*!>*}) < 

where M > is a bound for the sequence \\ipE n Wh 1 > an d such a bound exists because the sequence 
il>E n is weakly convergent in H (R). Then, by Rellich's compactness theorem, 

lim \\ipEn - iPe.\\l*({\x\<R}) =0 > 
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hence there exists uq such that 



\\ipE n ~ ipE. \\l*({\x\<R}) < TTfTTTTi , for n > n , 

2 \\ v IIl°°(m) 

and, consequently, 

\\VlpE n ~ VlpEt IIl 2 (R) < [|V|U«'(K)IIV'Sn ~ ^EM^dlxlKR}) + \\V\\l°°({\x\>R}) \\i>E n ~ ^E,\\l 2 (R) 

< - + — 2M < e, for n > n 
2 AM ~ ' ~ 

which proves VipE n ~^ V^Pe* ■ Now both terms in the bracket of the right hand side of (|3.7p converge 
in L 2 (M), and (— <9 2 + -En) 1 — ► (— <9 2 + E*)^ 1 in the space of bounded linear operators from L 2 (R) 
to H 2 (W), we obtain ■i/'E n convergent in H 2 (M). But _ff 2 (IR) is continuously imbedded in L q (M), for 
2 < q < oo, hence the limit of ipE n in -ff 2 should be the same as its limit in L g , 2 < q < oo. 
Therefore, (|3.6p implies that V'En converges to ipE* in .£f 2 (IR). Moreover, by passing to the limit in 
(l3~7l) . we get 

= {-d 2 x + E^-V^e^x) - a\^ E ,{x)\ 2p ^ Et l 

or equivalently, (ipE. t ,E*) is a solution of the stationary NLS equation (jl.3p . 

In addition, the eigenvalues of L + calculated at (ipE„,E n ) will converge to the eigenvalues of 
L + calculated at In particular, at (ipE«,E*), L + will have the lowest eigenvalue strictly 

negative, followed by as the second eigenvalue, see Lemma [2] and case (b) in Theorem [2j By 
Sturm-Liouville theory is a simple eigenvalue with corresponding eigenfunction being odd. Hence 
L + at (ipE„,E^) is invertible with continuous inverse as an operator restricted to even functions. 
The implicit function theorem applied to (II. 3p gives a unique C 1 in H 2 (R) branch of symmetric 
(even in x) solutions in a neighborhood of (ipE»,E*). By uniqueness, this branch must contain 
(ipE,E), Eg (Ei,E*). In particular: 

Jim \\tpE ~ ipE t \\m = 0, 
E-+E* 

and Lemma [3] is completely proven. □ 



We now return to the proof of Theorem O Based on (|3.3p and on the weak relative compactness 
of bounded sets in the Hilbert space i^ 1 (M) we can construct a sequence E n /* E* and ipE* £ 
such that (|3.4p and (|3.5[) are satisfied. It suffices now to show that at least for a subsequence (13. 6p 
is satisfied. Define 

= (3 ' 8> 

and the concentration function: 



<y-t 

see [U Section 1.7] for its main properties. Let 



ry+t 

t) = sup / \(p(x)\ 2 dx, (p G L" 
y m Jy-t 



u = lim liminf p(^f n ,t). 

t—>oo n— >-oo 
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Then, from (|3.3|) and part (i) of Theorem El the sequence \P n is bounded in and normalized 

in L 2 (M). According to the concentration compactness theory, see for example [5l Section 1.7], we 
have the following three possible cases. 

1. - Vanishing: if fj, = 0, then there is a subsequence ^ nfc convergent to zero in L 9 (M), 2 < q < oo, 
and, using (13.8|) and part(i): 

L q 

*pE nk — > 0, for all 2 < q < oo. 

Via Lemma [3] we now get 

rl>E nk ^ 

which contradicts the fact that liniE-^ HV'eJI^ = ^ 0, see part (i) of the theorem. 

2. - Splitting: if < pi < 1, then there is a subsequence the sequences v k , w k , z k G -ff 2 (R) 
all of them bounded in H , the sequence i/j 6M and the function v* G such that: 

fe t = UE n jL^n k =V k + W k + Z k , (3.9) 

supp Vfc n supp w k = $, lim dist(?/fe, supp u;&) = oo (3.10) 

Vk(--Vk) v*, v k (- - y k ) -)■ u*, 2 < q < oo, (3.11) 

|H|| 2 = /xiV* > 0, lim \\w k \\ 2 L2 = (1 - /i)iV* > (3.12) 

k— >oo 

z k ^ 0, z k 0, 2 < g < oo. (3.13) 

If, in addition, the sequence y k G M is unbounded, then, by possible passing to a subsequence, 
we have lim^—^oc, y k = ±oo. Both cases can be treated exactly the same so let us assume 
linifc^oo y k = oo. Fix a compactly supported smooth function 4> £ Co°0^)- Then, from the 
weak formulation (11.31) we have: 



= (<£'(•), (v k + w k + z k )'(- - y k )) + (</>(•), - + E nk )(v k + w k + z k )(- - y k )) 
+&((/>(■), \v k +w k + z k \ 2p (v k + w k + z k )(- - y k )). 

By passing to the limit when k — > oo, the terms containing v k will converge to v*, see (13. lip , 
the terms containing w k respectively z k will converge to zero due to (|3.10p respectively (|3.13p . 
and the term containing V(- — y k ) converges to zero due to lim^i^cx, V{x) = 0. Hence 

(</>',<) + + a((/),\v*\ 2p v*} = 0, 

and G ff 1 (M) is a weak solution of 

-v'l + E*v*+o\v*\ 2p v* = 0, (3.14) 

therefore a strong, exponentially decaying solution, with v'^ also exponentially decaying. Let 

t t \_ v *( x + Vk) - v*(~x + Vk) 
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£ k (x) is obviously odd in x, and the exponential decay of v* and its derivative implies: 
lim (v*(x + y k ),v*{-x + y k )) = 0, lim (v*(x + y k ),L + v*{-x + y k )) = 0, 

k— >oo fe— >oo 

where L + is calculated at i? nfc . Consequently: 
, lim ll&llfa = 1, 

fc— >oo 

lim (&., = -— lim (v*(x + y k ),L + v*(x + y k )), 

k— >oo Ul\* k—>oo 

where we also used f)3. 12|) for the first identity, and the fact that L + is invariant under 
the transformation x t-t —x because both ipE„ k (%) and V(x) are even in x, see also (jl.5p . 
Employing again the splitting (|3,9p , the convergence (|3,lip , f)3. 13|) , and the separation (|3.10p 
we have: 

lim (v*(x+yk),L + v*(x+y k )) = — v'l+E^+{2p+\)a\v^\ 2p v^) = 2pa l \v*\ 2p+2 (x)dx < 

fe^oo ' J R 

where we also used (|3.14p . All in all we have 



lim <£ fc ,L+£ fe ) = ^ / \v*\ 2p+2 (x)dx < 



fc— »oo 

which together with ^ k odd implies that for k sufficiently large the second eigenvalue of L+ 
at E n , must become strictly negative, in contradiction with Remark [3l 

If, on the other hand, the sequence y k E K is bounded, then, by possible passing to a 
subsequence, we have lim^oo y k = y* £ M. Let u = + y*). We will now rename u to be 
and we get from (|3.1ip 

v k — v u*, Vfc — >• «*, 2 < g < 00, (3.15) 

and (|3.12p remains valid. As before the convergence implies that v* satisfies weakly, then 
strongly: 

- v" + V(x)v* + E*v* + a\v*\ 2p v* = 0, (3.16) 
hence, both v*(x) and v^(x) decay exponentially in x. Consider 

jl = lim liminf p(w k ,t), where < Jl < (1 — p,)N*. 

t— >co fc— 5>oo 

If Jl > then w k splits: 

Wk = v k + w k + z k 

where the sequences v k , w k , z k E H 2 (W) have properties f|3. 10j) — (|3. 13|) . In particular there 
exist v* E H (K), ||£*||l 2 = A an d the sequence y k E M such that 

Vfc(- - Jffc) 5*, v k (- - y k ) 5*, 2 <q< 00. 

But now, the sequence is definitely unbounded since E supp and y^ is bounded, see 
(fBTTOl) . Using 

u*(x + y fe ) -^(-cc + yfc) 

Wi) - 
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we get as before 

lim <£ fe ,L + £ fc ) = ^ f \v*\ 2p+2 (x)dx < 

which contradicts the fact that the second eigenvalue of L+ should remain non-negative, see 
Remark [3l 

If, on the other hand, \x = then, by possibly passing to a subsequence, we have 

H 1 L q 

Wk — 1 0, Wk — > 0, 2 < g < oo. 
Moreover, from (EHJ) , (13TT31) and (IXT5I) we get 

i>E nk «*, ^£ nfc 2 < g < OO, 

which, via LemmaEl implies lim^oo HV'En — ^ll/f 2 = 0- The latter is in contradiction with 

\\v*\\l 2 = (iN* < iV* = lim \\ijj En \\ L 2, 

3. - Compactness If /i = 1 then there is a subsequence ipE n , a sequence ^ el and a function 
G H 1 ^), \\ip*\\ L 2 = iV* such that: 

^E nk (• - Vk) ip*, ^E nk (• - Vk) ->■ 2 < g < oo. 

But the sequence must be bounded, otherwise, by possibly passing to a subsequence, 
lim/^oo yt = ±oo. Both cases are treated in the same way, so let us assume lim^oo y k = oo. 
Then, by the symmetry of ipE n we have: 



fe t 0*0 +_fe t (~ x ) L\ A (x + y k ) + V* {-x + y k ) 
2 



and we get the contradiction: 



lim U Enk {x)\\i^ = lim H^llia, 

fe— >oo K Z fe— >oo 



where the last identity is a consequence of lim^oo (-0* (x + yk),ip*(—x + = 0. 

So, the sequence 6 R must be bounded, and, by possible passing to a subsequence, we 
have lim^oo y k = y*, and 

fe„ fe ipE, , ^E nk ~^ipE,, 2 < q < oo, 

where ipE*{x) = tp*(x + y*). Therefore the hypotheses of Lemma [3] are verified and part (ii) 
is proven. 

Theorem [3] is now completely proven. □ 

Remark 4. The proof of Theorem can be greatly simplified if one knows apriori that, with the 
exception of the two lowest eigenvalues, the spectrum of L + is bounded away from zero. 
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Indeed, in this case only the second eigenvalue of L + can approach zero, see Lemma [21 Hence L + 
restricted to even functions is invertible with uniformly bounded inverse. In particular 

9eiPe = -L^ipE, \\8eiPeWl 2 < KW^eWl 2 

The argument in Theorem [2] can now be repeated to show directly that ipE has a limit in L q , q > 2, 
and the proof is finished by applying Lemma 

Corollary 1. Under the assumptions of Theorem^ we have 

(i) the second eigenvalue of L + (ipE, E), E < E*, denote it by X(E), and only the second eigen- 
value approaches as E /* E* ; 

(ii) if in addition p > 1/2 and the derivative of the second eigenvalue of L + satisfies: 

then the set of real valued solutions of f| 1 . 3 1) in a H 2 x IR neighborhood of (ipE*, E*) consists 
of exactly two curves of class at least C' 2 ^ -1 intersecting only at (ipE t ,E*). 

Proof. For part (i), Theorem [3] part (ii) implies that L + depends continuously on E G [Eq,E*], 
hence its isolated eigenvalues will depend continuously on E. In particular, at E*, L + will have the 
lowest eigenvalue strictly negative, followed by as the second eigenvalue, see Lemma [2] and case 
(b) in Theorem [2 If any other continuous branch of eigenvalues of L + approaches zero as E /* E*, 
then becomes a multiple eigenvalue for L + at E* in contradiction with Sturm-Liouville theory. 

Part (ii) is a standard bifurcation result which uses Lyapunov- Schmidt decomposition and Morse 
Lemma, see for example [24J. More precisely, we continue working with the functional F((j),E) : 
H 2 (R) x R i->- L 2 (R) given by (|2.ip which has the Frechet derivative: 

D^F((j), E) = -8 2 x + V{x) + E + (2p + l)o-\<J)\ 2p = L + {4>, E). 

For p > the functional is C 1 while for p > \ it is C 2 . By part (i) D^F(tpE,,E^) has zero as a 
simple eigenvalue. Let 0* be the L 2 -normalized eigenfunction of = D^F^ipE^, E*) corresponding 
to its zero eigenvalue. Then is a Fredholm operator of index zero with KerL* = span{0*} and 
RanL* = [KerL*] -1 . Let 

P\\<t> = (0*^)^20*) p±4> = <P - (<P*,<P)l 2 4>* 

be the two orthogonal projections associated to the decomposition L 2 = KerL* © RanL*. By the 
standard Lyapunov- Schmidt procedure, see for example [24], we get the following result. 

Lemma 4. (Lyapunov-Schmidt decomposition) There exists a neighborhood W C H 2 x R of 
{ipE t i E*), a neighborhood U Cl 2 of (0, E*), and an unique C l map h : U i— > L 2 n {4>*}~ L such that 
for any solution (4>,E) £ W of F((j),E) = there exists a unique a £ M. satisfying: 

(a,E)eU, 4> = ip Et + a4>* + h{a, E) 
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and 



F {l (a,E) = {fa,F(il> E , +a<p* + h(a,E),E)) = 
In addition, for all (a, E) € U we have: 



h(0,E*) 

da- {a > E) 
dh 



dE 



(a,E) 



-(P L L+)- 1 P L L+^ J 
-(P ± L+)- 1 P ± [V^ +h], 



where L + = L + (ip Et + acfi* + h(a, E),E). 

Moreover, if p > \ then h is C 2 and for all (a, E) G U we have: 



d 2 h 
dadE 
d 2 h 



dE 2 



(a,E) 
(a,E) 



-(P ± L+)- 1 P ± (a a L 
dh 



dh 
da 



-(P ± L + )- 1 P 1 
-(P ± L + )- 1 P 1 



m r \^h 

d~a + {9aL+) dE. 

dh dh 
9E + (9 E L + )— 



where 



(dh \ 

dh dh ( dh 

(d E L + )— = — + (2p + l)2pa\^ Et + afa + /f^signfe + afa + h) ( — 



(3.17) 

(3.18) 
(3.19) 

(3.20) 



(3.21) 
(3.22) 
(3.23) 



So, the solutions of (|1 .3j) in W are given by the solutions of ()3.17p . From (|3.18|) and F{ip Et , E*) 
we have i*li(0, E*) = 0. Differentiating (I3.17D we get: 



dE 



dh 



dF \\ 

~dE 



(a,E) 



i , i t dh . . , dh. 

m WE. + a-9* + h + L +^) = ° + \<P*> L +~q^i 



In particular, from (|3.19|) - (|3.20|) . L*^* = 0, and L* self adjoint, we have 

dFu 



da 
dF\\ 

~dE 



■(Q,E*) 
(0,E*) 



0, 
0. 



(3.24) 
(3.25) 

(3.26) 
(3.27) 



Since the gradient of Fn at (0, E*) is zero, the number of solutions of (|3. lTj) in a small neighbor- 
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hood of (0, E*) is determined by the Hessian at (0, E*). For p > i we can calculate: 

d 2 F\\ dh d 2 h 

-^J-(a,E) = {fa, {2p+l)2pa\^ Ef +afa + h\ 2p - l s\gn^ E ,+afa + h){fa + —) 2 + L + ^), 

d 2 F\\ dh dh 

■^-^(a,E) = 1 + (fa, (2p+l)2pa\i} Et + a^* + fr| 2p_1 sign(^ + afa + h)(fa + —)—) 

+(< ^' L+ ^dE ] 

d 2 Fu _ ,, dh , ,,o. i . ,. . w ah„ „ <9 2 /i, 



(a,£) = ^ + (2p+l)2pa|fe +a0* + /i| 2p " 1 sign(^, +a^ + /i)(—) 2 + L H 



In particular 

= (0., (2p+l)2paVg- 1 ^) + (^, L*0> =0, (3.28) 

d 2 F\\ Bh d 2 h 

;(0,E*) = l + (fa,(2p + l)2pa^- 1 fa—} + (fa,L,- 



dadE y ' *' vr " v ^ / ^ ^aE' \™ * dadE 1 

= l + (2p + l)2pa [ ^tf^-foEjdx (3.29) 
Jr aE 

d 2 F\\ dh dh d 2 h 

^j(0,£*) = (0*. ^> + <<^ (2p+l)2pa^%- 1 (—) 2 } + (fa, L„— )=0, (3.30) 

where we used (|3.18[) - (|3.20|) . the fact that is even while fa is odd, the fact that h and its partial 
derivatives are in {fa}' L , and the fact that L* is self adjoint with L*fa = 0. We next show: 

f dh d 2 F 

A'(£*) = 1 + (2p + l)2pa jf ^P-^l—tO, E*)dx = q^(0, E.). (3.31) 

The second eigenvalue of L + (ipE, E) along the C 1 symmetric branch (ipE,E) satisfies the equation 

L+{i/> E , E)4> E = X(E)faE, Ue\\v = 1- 
Differentiating with respect to E, we get: 

4> E + (2p + Wp*^- 1 ^ + = AWE + A(^)^, 

and, by taking the scalar product with <f>E, we obtain 

1 + (2p + l)2pa / ^T^l^dx = X'(E). (3.32) 

Using the continuous dependence of the spectral decomposition of L + (ipE, E) with respect to E G 
(Eo,E*], we have lime ^e, \\</>e — ^Wh 2 = 0- Moreover, from (|2.5p and (|3.20p we have 

E% " = Jim \\L-^e - L^ Ef \\ H2 = 0. 

Passing now to the limit E /* E* in the identity above we get (|3.3ip . 
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Prom GEggD-flESJ we have 



A' (22*) 
A' (22*) 



Since by hypothesis A' (22*) ^ 0, the Hessian of 2*y is nonsingular and negative definite at (0,22*), 
hence by Morse Lemma, see \24\ Theorem 3.1.1 and Corollary 3.1.2], the set of solutions of 
Fu(a, E) = in a neighborhood of (0,22*) consists of exactly two curves of class C^]" 1 inter- 
secting only at (0, 22*). This finishes the proof of the corollary. □ 

Unfortunately, for 1/2 < p < 3/2, the corollary does not guarantee that the two curves of 
solutions are C 2 which turns out to be necessary for determining their orbital stability with respect 
to the full dynamical system (II. lh . However a more careful analysis of equation (13.171) recovers the 
full regularity and stability of these solutions: 

Theorem 4. Let a < 0, p > 1/2, and consider the symmetric (even in x) branch of solutions 
(ipE,E) of (|1.3p which bifurcates from the lowest eigenvalue —Eq of L$. Let (22o,22*) be the maximal 
interval on which this branch can be uniquely continued. Assume 22* < oo, and 

A'(22*) = lim ^(22)^0, (3.33) 
E dE 

where A(22) is the second eigenvalue of L + (ipE, 22). Then, the set of real valued solutions (<j),E) G 
H 2 xl of the stationary NLS equation (jl.3p in a small neighborhood of (tpE* , 22*) G H 2 x R consists 
of exactly two C 2 curves intersecting only at (t/^,22*) : 

(i) the first curve can be parameterized by 22 h-> (ft = tpE, 22 G (22o, 22* + e) for some small e > 0, 
it is a C 2 continuation past the bifurcation point E = E* of the symmetric branch, it has ipE 
even for all 22 and orbitally unstable for 22 > 22*; 

(ii) the second curve is of the form (<p(a), 22(a)), a£i small, where the parameter can be chosen 
to be the projection of (ft — ipE, onto kerL* = span{(/>*} i.e. Be > such that for \a\ < e : 

E = 22(a) = 22* + ^a 2 + o(a 2 ), (ft{a) = ip E , + a0* + h(a), (3.34) 

where h(a) = 0(a 2 ) G {(ft*}' 1 , and 

q = - 2p % E ^ a2 (^<^r 2 ^ 2 - 2p(2 P+ i)^rv* 2 ,^-vrv* 2 >L 2 ) (3.35) 

with L* = 22*), along this curve (ft is neither even nor odd with respect to x, and is 

orbitally stable if 

Q > 0, and R > 



and orbitally unstable if 
where 



Q < 0, or Q > and R < 0, 



R = e™e~-^ = 2 ^^ + N ' {E * ] > N ( E ) = WM\h- (3-36) 
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Proof. We continue to rely on the Lyapunov-Schmidt decomposition described Lemma U but we 
remark that ipE* = linig /-e^'Pe given by Theorem [3] is even, since ipE, E < E*, were even. Hence 
the Prechet derivative of F(ip,E) with respect to ip at (t/>_e„ , -E 1 *) : 

D^F^e^E.M =L*[cj>} = (-d 2 x + V + E*)<j) + (2 P + l)\iPE t \ 2p (p, 

transforms even functions into even functions. By Corollary [T] part (i), is the second eigenvalue 
of L* hence its eigenfunction 0* is odd via Sturm-Liouville theory. Consequently 

L* | even '■ H e ven I— >■ L &ven 

is an isomorphism. Implicit function theorem for F(ip,E) = at (ipE*,E*) implies that the set 
of even, real valued solutions of (|1,3|) in a neighborhood of (i()e*,E*) consists of a unique C 1 
curve, E \-t tpE- Moreover i\>e > 0, because tpE is also the eigenvector corresponding to the lowest 
eigenvalue of 

L_ = -d 2 x + V + E* + \iiEj 2p , 

and, consequently, F(E) = F(x/je,E) becomes C 2 even for p = 1/2. Hence L + = D^F is C 1 in E 
along this curve. Differentiating once F(ipE-,E) = we get: 

d4>E ,j s.-i , 

and the curve E i— > tpE is C 2 because the right hand side is C . Moreover, the second eigenvalue of 
L + along the curve, X(E), is C 1 in E, and, because: 

A(£*) = 0, ~^( E *) < °. 

see Corollary [J part (i) and (j3.33|) . we deduce that X(E) < 0, E > E*. 

For E > E* we now have that L + has two strictly negative eigenvalues while L_ has none along 
the symmetric (even in x) branch of solutions of (ll.3p . These imply that exp(iEt)ijjE is an orbitally 
unstable solution of (jl.ip . see [8], and finishes part (i) of the theorem. 

For part (ii) we will rely on the curve of solutions discovered in part (i) to do a sharper analysis 
of equation (13.17j) compared to the one provided by Morse Lemma. The solutions (ipE, E) of (I1.3P 
discovered in part (i) satisfy a = ((/>*, ipE) = because ^* is odd while ipE is even. Consequently 
one set of solutions of (I3.17P is given by a = 0, in particular: 

ipE = ipE*+h(0,E) >0 (3.37) 
F ]{ (0,E) = {fa, F(tp E *+h(0,E),E)) =0 (3.38) 

Hence, "a" can be factored out in the left hand side of equation (I3.17j) . and, solutions a ^ of this 
equation satisfy: 

g(a,E)=0 (3.39) 



where 



We will show that: 



r F||(a,E)-F||(0,E) 

g(a,E) = { aFi a , na^u , . 

I ifa = 
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(a) g(0,E*) = 0, 

(b) g is C 1 in a R x R neighborhood of (a = 0, E = E*), 

(c) §a.(0,E*) = \>(E*)^Q, §|(0,^) = 

Implicit function theorem will then imply that the set of solutions of (|3.39j) in a neighborhood of 
(0,E*) consist of a unique C 1 curve: a i-> -E(a), |a| < e for some e > 0, with E(0) = E* and 

^ (0) = -|(^l = o. 

■fa" f$(0,fi.) 
Then we will show that the curve a \— y E(a) is C 2 with 

0(0) = g (3.4i) 

where Q is given by (13.35H . Hence the nonsymmetric solutions of (11. 3p . are given by: a i— > 
(<f>(a),E(a)) with: 

E = E(a) = E* + ®a 2 + o(a 2 ), 0(a) = + a<£* + %, £7(a)) > (3.42) 

and a i— )■ </>(a) is C 2 (from R with values in ii 2 ) because from (|3.42p and f)3. 18|) : 

^p-(a) = ^ + ^+ E'{a)^ = 0* - (P ± L + y 1 P ± [L + ^ + E'{a){4>* + h)] 
da da dE 

and the right hand side is C 1 in a because E, h are C , and L + is also C 1 when calculated along 
(0(a) > 0, E(a)). The orbital stability of this curve of solutions follows from the theory developed 
in |B], [2] and: 

Ai(o) = -\'(E*)Qa 2 + o{a 2 ) (3.43) 
||</»(a)|| 2 2 = N(E*) + ^(2\'(E*) + QN'(E*))a 2 + o(a 2 ) (3.44) 

where Ai(a) is the second eigenvalue of L + (cf>(a), E(a)). Indeed, for Q < 0, ()3.43p shows that L + has 
two strictly negative eigenvalues while we know that L_ has none, hence the result in [8] implies 
(j>(a) is orbitally unstable. For Q > 0, L + has exactly one strictly negative eigenvalue, L_ has none, 

and the result in [9] implies that 4>(a) is orbitally stable if > and unstable if < 

which via (|3.44j) and (|3,42p is equivalent to it! > 0, respectively R < 0, where it! is defined by (|3.36p . 

Now AMU) follows from flS32D , (1331]) and properties (fBTTH^ - (f3T23]) of /i(a, E), while QgaSD follows 
from the same relations by differentiating twice the equation for the second eigenvalue: 

L + (^ E (a),E(a))[<p 1 (a)} = Ai(a)0i(a), ||<Ma)|| L 2 = 1. 

It remains to prove (a)-(c) and ()3.4ip . (a) follows from ()3.40p and ()3.26p . For a / (b) follows 
from F being C . For a = it suffices to prove: 

dg dg d 2 Fu 

lim —(a,E) exists, and hm —(a,E) = (0,E ). 
a^o da ajto oE oEoa 
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Note that g E Q a (0, Eq) exists and it is continuous in Eq because the partial derivative with respect to 
E of the right hand side of (|3.24p is the derivative along the symmetric branch ip E = tpE, +h(0, E) > 
which we already know is C 2 . We will prove that the first limit exists, the argument for the second 
is similar. In what follows we use the shortened notation: 



L+(a, E) = L + (ip Et + + h(a, E),E) = D^F{^ E , + a<t>* + h(a, E),E). 



For o^Owe have 

J-(a,E) = Afa, F(i> Et +a<f>*+h(a,E),E)-F{i> Et +h(0,E),E))+~{<f>*, L + {a,E)[<f>*+—(a,E)]). 

oa a z a oa 

We add and subtract ±(</>*, L+(0, + §|(0,£)]) to get 

lim Tr(a,E) 
a^o da 

a^0,E^E 



Inn 

o 

a-+0,E^>E 



F\\ (a, E) - F,| (0, E) + a(^, L+ (0, E) [</>* + § (0, E)}) 



a^o a 2 



+ Um (<k, L + (a,E)[^ + f a (a,E)]-L + (0,E)[^ + ^(0,E)}) = ^ + 

a^o a 
a->0,E->-Eo 

The differential form of the intermediate value theorem gives: 

h = _i lim (0*, MgVgH^ + j(gVg)] - MMQ[^ + jgVg)]) 

2 a' 

where a' = a' (a, is between and a. So Ii = —5/2 provided I2 exists. 

For the limit I2, when p > 1/2 we can use F is C 2 hence L + is C 1 and h is C 2 to get: 

dh d 2 h 
h = d a L+(0,E )[^ + —(0,Eo)] + L + (0,E )-^(0,E )) = 0, 

see ()3.2ip for an explicit expression for both d a L + and and notice that they are both even at 
(0, Eq) while is odd. When p = 1/2 we rewrite I2 as two limits: 

/2 = Hm (0*. (L + (a,E) - L + (0,E))[^ + f («,£)]) 
a#o a 

a->0,.E->-Eo 

+ hm — — = J3 + J4 

a#o a 
a-)-0,E-^E 

But 

(L+(a,E)-L + (0,£;))H =MW>E. + a<f>* + h(a,E)\ -\^ Em + h(0,E)\)v 
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and by —\a — b\ < \a\ — \b\ < \a — b\ and for each x G 



\h(a,E) - h(0,E)\(x) < 



dh , . _,. . , 
— (a',E)(x) 



\a\, for some la'l < lal 



we get that the integrand in the expression for the limit ^3 is bounded by: 



Mx)(L+(a, E) - L+(0, E))[cf>* + f (a, E)](x) 



< 



(x)+6J(x) 



d-a {a > E) 



where the right hand side can further be bounded by an integrable function, since h(a, E) G H 2 is 
C 1 in a neighborhood of (0, E*). 

In addition, since (a,E) h-> h(a,E) G H 2 (R) ^ L°°(M)nC(E) is continuous and ip E * +h(0, E )) > 
0, see (|3.37p . we get that for each i£l there exists a 5(x) > such that ipE* + &4>* + h(a, E)) > 0, 
for \a\, \E — Eq\ < 5. So, we have the pointwise convergence: 



hm 

a^0,E^E 



\ipE,+a^+h(a,E)\(x)-\ifj E ,+h{0,E)\(x) 



= ^ oM*) + (h(a,E) - h(0, m x) = + dh x ^ e 

a^o a da 

a^0,B->E 

which combined with the Lebesque Dominated Convergence Theorem implies 

dh 

h = (</>*, 2a(^ + 7 -(0, J E )) 2 ) =0, 

since 4>* is odd and (0* + ^{0, Eq)) 2 is even. 

Similarly, from (|3.19p we get the pointwise convergence: 



f(q,£)|(x)-f(o,£)Qr) 



-(P ± L+)- 1 P ± (a a L)[^ + — (0,^o)](^)Vx G K, 

aa 



and, again by Lebesque Dominated Convergence Theorem: 



dh 



h = ~M<t>*, L + {^E G ){P L L + )- 1 P 1 _[4 )if + — (0,E )} 2 ) = 0, 

since <f)* is odd while the other factor is even. 
In conclusion, for p > 1/2 we have 

da 1 dh d 2 h 

lim -£(a,E) = -{cj ) „ d a L + (0,E )[<j)* + —(0,Eo)]+L + (0,E )-^ I (0,E ))=0. 

a^,0,E->E 



A similar argument shows 



dg 



d 2 E 



a J m dE^ E) dEda 
a^0,E^E 
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where 

d 2 F\\ Bh dh d 2 h 

^±(0,E ) = l+(2p+l)2jw7(^, (^E t +h(0,E )) 2p ~\^ + — (0,E )) — (0,E )+L + ^^(0,E )). 

So, g(a,E) is C 1 with 

g(0,E.) = 0, and ^(0,E,) = 1 + (2p + l)2pa J ^^^(Q, E,)dx = \'(E,) ^ 

see (|3.3ip and (|3.33p . Hence the set of solutions of g(a, E) = in a neighborhood of (0, E*) consists 
of a unique C 1 curve a i— > -E(a) with .£7(0) = E* and 

-(0) = -f-(0,^) = 0. 

The curve is in fact C 2 , because for a ^ 0, p > 1/2 we have that g(a,E) is C 2 while for a / 
0, p = 1/2 a similar argument as the one above involving pointwise convergence to i^E( a ) = 
ipE* + + h(a, E(a)) > can be employed. For a = we have from definition of the derivative: 



E"(0) = lim - f (a '^ (a)) = —1— lim ^(a,£(a)), 



where 



lim ^(a,E(a)) = \{2p + l)2p{2p - l)^ 2 , V^T^*) - ^p ^ \)\2pf ^K^V €) ■ 

The last identity follows from the same argument as for lim^o, §f ( a >^) with the only differ- 
ence that for 1/2 < p < 1 Lebesque Dominated Convergence Theorem requires: 



\il>E. + a0* + E(a)) \ 2p ~ 2 \(f>* \ £ L 2 f| L°° . 



But since both 0* and iPe{o) = V's* + o</>* + h(a,E(a)) > are solutions of the uniform elliptic 
equations: 

= 0, L-ipE(a) = 0, 
we have via upper and lower bounds for uniform elliptic equations: 

10*0*01 - C(5)exp(-y / E^ - S\x\), for 8 > 0, and^(a) > C(e) exp(- y^* + e|x|), for e > E(a)-E*. 
Hence, for 1/2 < p < 1 we get: 

|^(«)| 2p ~ 2 <Cexp(-e|x|) 

for any e > such that: 

E*+e< 



2-2p 

This finishes (a)-(c) and (|3.4ip . consequently Theorem [H □ 
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We note that Theorem 2] shows that the branch of symmetric states goes through a pitchfork 
bifurcation at (iJje*, E*), which is classified, based on stability analysis, as supercritical when Q > 
and R > 0, and subcritical when either Q < or Q > and R < 0. We do not have a general 
result determining which one occurs except for the double-well potentials with large separation, e.g. 
(jl.4p with s large: 



Corollary 2. Consider the equation (I1.3P with a < and potential of the form: 

V = V,(x) = V {x + s) + Vq{-x + s), x E R, s > 0. 

Assume Vq satisfies (HI), (H2) and (H4). T/ien i/iere exists > suc/i i/tai /or a// s > s* the 
branch of real valued, even in x solutions (tpE,E) bifurcating from the lowest eigenvalue of —A + V s 
undergoes a pitchfork bifurcation at E = E*, where lim s _>. 00 E* = Eq. Moreover, the asymmetric 



branch emerging at the bifurcation point is orbitally stable if p < p* 



3+V13 



and orbitally unstable 



if p > p*, while the symmetric branch (i/je,E) continues past the bifurcation point but becomes 
orbitally unstable. 

Remark 5. We remark that the case p = 1 has already been obtained in |i6] /. In what follows 
we present a more direct argument to obtain the same result for all p > 1/2. The argument can be 
easily adapted to higher space dimensions, see [20], and to the case a > in which we can rigorously 
show that a pitchfork bifurcation occurs along the first excited state ( the branch bifurcating from the 
second lowest eigenvalue of —A + V s ). This result has been predicted in [30]. 



Proof. Under hypotheses (HI), (H2) and (H4) for Vq, it is known that the spectrum S of —8% + V s 
has the lowest eigenvalues —Eq^ s < —E\ tS satisfying 

lim \E ks -E \ =0, k = 0,l, and 3 s*, d > : dist(£o 8 , S \ {E x „}) > d Vs > s*. (3.45) 

Moreover, the normalized eigenfunctions iposj^is corresponding to these eigenvalues satisfy: 

ip (x + s) + ip (-x + s) 



lim 

s— >oo 



lim 

s— >oo 



V2 

ipo(x + s)- ip {-x + s) 



V2 



H 2 



H 2 



0, 
0, 



(3.46) 
(3.47) 



where ipo is the eigenfunction of —8% + Vq corresponding to its lowest eigenvalue —Eq, see [TBI 
Appendix] and references therein. Proposition [1] shows that for each s > 0, a unique curve (if)E,E) 
of real valued, nontrivial solutions of 



-dl4> + V s 4> + E<f> + a\ 



|2P, 



0, 



bifurcates from (0, Eq) and can be parametrized by a = (tPo^^e) , i-e. there exists e > 0, such that 
for Id < e we have: 



E 

ipE 



^o-o-||Vo, s ||^ + 2 2 |a| 2p 



+ 0(|a| 4p ), i.e. \E-E -a 
aVo,s + 0(|a| 2p+1 ), i.e. \\^ E - a^ , s \\ H 2 < C 2 \a\ 2p+1 . 



l ||2p+2 I 1 2p I 



< Ci|a| 4p (3.48) 
(3.49) 
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Moreover, i/je( x ) is even in x. We will rely on the fact that the distance between the two lowest 
eigenvalues of d 2 + V s converges to zero as s — > oo, see (|3.45p . to show that there exists s* > such 
that for each s > s* the second eigenvalue A(a, s) of L + (ipE,E) must cross zero at a = a* = a*(s), 
where < a* < s, and 

lim a* (a) = 0. (3.50) 

Moreover, for each s > we have 

< "P<0. (3.51) 

Hence the hypotheses of Theorem [J] are satisfied and a pitchfork bifurcation occurs at E* = E{a*) 
for each s > s*. We will also calculate Q = Q(s), R = R(s) at these bifurcation points and show 
they are continuous on the interval s £ [s*,oo) with 

lim a 2 -^{s)Q{s) = -a^(2p + l)(p + l)||Vo|lS+ 2 2 (3.52) 

lim a^- 2 ( s )fl( s ) = 2P(-p 2 + 3p+l) 

s — -a(2p+l)(p + l)p||^o||?i 2 2 

Therefore, by choosing a larger if necessary, we have Q(s) > 0, for all s > s* and, if p < p* = 
3+ ^ then R(s) > 0, for all s > while if p > then _R(s) < 0, for all s > s*. The proof of the 
corollary is now finished. 

It remains to prove (|3.5U|) - (|3.53|) . They follow from rather tedious calculations involving the 
spectral properties (|3.45p - (|3.47p . and bifurcation estimates (|3.48p - (|3.49|) . We include them for 
completeness. First we note that there exists s* > such that the estimates (|3.48|> — (|3.49j) are 
uniform in s > s*, i.e. the constants C±, C2 can be chosen independent of s > s*. The reason 
is that the estimates rely on contraction principle applied to the operator K ((f), E, s) : H 2 1— > H 2 
given by: 

K{J>, E, s) = I - (-8 2 x + V S + E^ s )- l P^ s {-d 2 x + V s + E + (2p+ l)aH 2p ), 

where P± s denotes the orthogonal (in L 2 ) projection onto ^oV Since the spectrum of (— d 2 + 
V s + Eo iS )P± tS restricted to even functions remains bounded away from zero for s sufficiently large 
and V s : H 2 1— > L 2 is uniformly bounded, we can choose the Lipschitz constant for K, hence the 
constants Ci, C2 above, independent of s. 

Now, if X(a,s) denotes the second eigenvalue of L + (ipE, E), then we have: 

X(a, s) = E 0)S - E 1)S + ^(0, s)(E - E 0>s ) + 0(\E - E^ s \ 2 ), (3.54) 

where, as before, the constant hidden in the 0(\E — Eo tS \ 2 ) term can be chosen independent of s 
for large s. Using now (|3.32p we get: 



dX 

— (0,s) = hm 
dE a^0 



l + (2p+lV / 0|A|^|2p d J =1 _ 2 P+ 1 r tf (x)1%,(x)dx 

Jr dtj \ \m, s \\L2 P +2 



where we relied on the expansions (|3.48p . ()3.49p . and on the continuous dependence with respect 
to a of the spectral decomposition of L + which implies that the eigenfunction <pE corresponding 
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to the second eigenvalue of L + converges to ipi s . Moreover, using the fact that the expansions of 
ipE, <t>E are uniform in s for large s, we get from (|3,46p - (|3.47p : 

lim a~ q [ ^%(x)(f) E k (x)dx = 2 1 - k - q / 2 \\^ \\ q I tl k 2k , for all q > 0, and fc = 1,2, ... , (3.55) 



a->0 

which in particular implies 



dX 



(3.56) 



lim -r-(a, s) = —2p < 0. 

s-i-oo dE 

a— >o 

From (|3.45p . (|3.54p and (|3.56p we now have 

lim X(a,s) = + 2pa2-P\\M%t+M 2p + 0{\a\^). 
Hence, because a < 0, there exists s* > and e > such that: 

A(e, s) < 0, and s) < —p < 0, for all s > s* and \a\ < e, a / 0. 

clE 

Using now A(0, s) = Eq jS — E\^ s > 0, s > we infer that for each s > s*, A(a, s) changes sign 
exactly once in the interval a G [0, e] at 



a = a*(s) 



-5b,s — Ely 
-A'(0,s) 



Hence ([330]) and (I3~5TD hold. 

We now compute Q and i? in the limit s — > oo, relying on (|3.50p . In formula (13.35P which defines 
Q we have already showed: 



lim \'{E*) = lim -^-(a,s) = —2p, 

a— >0 

and 

lim (^,Vt" 2 ^) = lim a 2 - 2 ? / ^(^^(xjdx = 2"' ||Vt>|g£? a 
see (I3.56P and (I3.55p . For the remaining scalar product we use 

= L + (ipE*,E*)TpE t = cr2pV 2 £ +1 
which for even functions is equivalent to: 

-1„/2p+1 



(3.57) 



* 2pa^ * 



We note that from (|3.55p we have 



tit 1 ?. = 



+ 



4k 



>E.\\h 



« +1 . A 



E.\\ L 2 S 
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hence 



/ -V-l 2 r -l 2- ' - '"/'-J r-l -'/'+!> i ^ /r-l/ 2 P+ 1 /2p-l / , 2 ^B, 



Passing to the limit when s — > oo and using (13.55h we have: 



Because the first two and only the first two eigenvalues of approach zero as s — > oo we need to 
expand the quadratic form in the last limit in terms of the associated spectral projections. Since 
the quadratic form involves only even functions and the eigenfunction corresponding to the second 
eigenvalue of is odd we only need to worry about the projection onto the eigenfunction £* 
corresponding to the first eigenvalue \q(E*) of L*. We have: 



The latter can be calculated via L'Hospital where, as in (13,561) . 

lim Ao(-E*) = lim ^^-(a,s) = —2p < 

a— >0 

and the derivative of the denominator converges to zero. All in all we get (|3.52p . 
Finally, to compute lim^oo R we use the definition (|3.36|) . We have 

/dtpE ( dE\ _1 , x _ 2a 1 +_0(\a*\ 2 P +1 ) 

,2p+2 

where we used (|3.48p - (|3.49|) . Consequently 

2 p 



d^ E (dE\~ l 2a,+0(K| 2 P+ 1 ) 



lim al p - 2 N'(E*) 



"HIV'O 11^+2 



which combined with (|336]l . (l337j) and (13321) gives (|333j) . □ 
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4 Behavior of the symmetric and asymmetric states for large E 



In this section we show that if the branch of symmetric states (i^e,E) bifurcating from (Q,Eq) can 
be uniquely continued on the interval E G (Eq,oo), i.e. case (a) in Theorem [2] holds, then, modulo 
re-scaling, this branch must bifurcate from a nontrivial, even solution of the constant-coefficient 
NLS equation: 

-u'{ y0 (x) + o\u 00 \ 2p u 0Q (x) + u 00 (x) = Q, UooeH 1 . (4.1) 
Since the above equation has exactly one such solution: 

/ 1 + p\ i 
Uoo{x) = ( J sechf(px), (4-2) 

we infer essential properties of the branch (ipE,E) via bifurcation theory. In particular we show 
that when V(x) has a non-degenerate local maximum at x = then L + computed at (ipE,E) has 
two negative eigenvalues for E large, contradicting Remark [3j This finishes the proof of our main 
theorem. 

However the arguments developed in this section tell much more about all solutions (i/je,E) 
of the stationary NLS equation (|1.3p for large E. Certain scaling of the ||V , s||z,2, ||V>_e||l2 P +2 and 
||VtAe||l 2 norms f° r large E emerges from Theorem Combined with the concentration compact- 
ness arguments, these norms imply that the stationary solutions either bifurcate from solutions 
(14. 2p translated to be centered at a critical point of V(x), or centered at infinity, see Remarks [7] 
and [H The stationary solutions bifurcating from a finite translation of (|4.2p are localized near a 
critical point of V(x), and if the latter is non-degenerate the orbital stability of these solutions can 
be determined, see Theorem [6] and Remark [6l 

Theorem 5. Let a < 0, and consider a C 1 branch of stationary solutions (tpE,E) for E £ (£"1,00). 
IfV(x) satisfies 

V{x) 6 L°°(R), lim V(x) = 0, xV'{x) E L°°(R), (4.3) 

\x\— >co 

then 



(i) there exists < b < 00 such that 



II / ii 2 p+ 2 

lim ^%±i = b (4.4) 

_ -_ll±l b (4 . 5) 

E->co E -~ 5 2 p+l 

lim " V t'! i2 = ^^6, (4.6) 



and, after the change of variables: 



u E (x) = R l '^ E {Rx), R = ^ E + v{Q y ( 4 - 7 ) 
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ue satisfies: 

lim = —ob > (4.8) 

E— >oo 

-u" E {x) + R 2 (V(Rx) - V(0))u E (x) + u E (x) + a\u E \ 2p u E (x) = 0. (4.9) 

(ii) if in addition ^Pe(x) is even in x and L + computed at (ip E ,E) has exactly one negative 
eigenvalue for all E £ (E\,oo) then u E defined above converges: 

lim \\u E — itoollif 2 = (4-10) 

E— >oo 

where Uoq ^ satisfies (j4. 1 j) and is given by (|4.2p . 

Before we prove the theorem let us note that it implies that, for the case (a) of Theorem [2j 
the branch of symmetric states, under the re-scaling (14. 7h : (ipE,E) i— > (ue,R), bifurcates from the 
solution (ue = UoojR = 0) of equation (|4.9p . where Uqo is given by (|4.2p . This bifurcation can be 
analyzed in detail: 

Theorem 6. Consider xq G R, £/ie equation: 

G{u, R) = -u" + R 2 {V(Rx + x ) - F(x ))u + n + a\u\ 2p u = 0, (4.11) 
and the solution (u = Uoo,R = 0), with given by (14.2p . 

// V is differ entiable at xq but V'(xo) ^ then the set of solutions of G(u,R) = in a 
neighborhood of (itoo>0) £ -H" 1 x 1R zs given by R = and translations of : 

{(«»(■ -«),o) : set). 

(m,) IfV is twice differentiable at xq and xq is a non- degenerate critical point ofV, i.e. V'(xq) = 
0, V"(xq) 0, then the set of solutions of G(u, R) = in a neighborhood of (n^, 0) G -ff 1 x R 
consists of two orthogonal C 1 curves: 

{(uoo(- - s),0) : s6l) and : i?6l, |22| small}, 

where 

hnWh = \M\h + - I) ^ / (x )||xn oo (x)||| 2J R 4 + O ( J R 4 ). 

Moreover, if xq is a local maximum ( respectively local minimum ) for V then L + computed at 
(ufj,R) has exactly two (respectively exactly one) negative eigenvalues. 

We now outline a few remarks. 

Remark 6. The last two theorems combined with Remark [3] finish the proof of our main theorem. 
Moreover, for the generic potential (ll.4p with s > s*, which has exactly two non- degenerate local 
minima at x = ±xq and one non- degenerate local maximum at x = 0, the above theorem gives two 
more branches of asymmetric states: 

E = R- 2 - V(±x ), ^e(x) = R-pur (^P-) (4.12) 
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localized near the two minima. Both branches are orbitally stable for p < 2 and orbitally unstable 
for p > 2. This is because the operator L + (^e,E) has exactly one negative eigenvalue, see part (ii) 
of Theorem [SJ and, according to the general theory in the sign of OeW^eW^ determines the 
orbital stability. Using (|4.12p we have 

d E UE\\h = (E + V{±x ))^\\u R \\ 2 L2 -\{E + V{±x Q )) 1 v- 2 d R \\u R \\ 2 L2 



1 3 
Ev 2 



Hence ^ll^jsllxa > 0, for p < 2, implying stability, OeW^eW^ < 0, for p > 2, implying instability, 
while for p = 2 we have 

Ue\\ 2 L 2 = WMh = hooWh -\v"{±x Q )\\xu 00 {x)f L 2R /i + o(R% 
i.e. \\iPeWl2 is increasing with E ~ R~ 2 , implying stability. 

Remark 7. The proof of Theorem [5] part (ii) shows that in the absence of hypothesis tpE even, the 
re-scaled ue may concentrate at ±oo, i.e. 

3yE G M such that lim yE = ±oo and lim \\ue(- — Ue) — ^oo||j3"2 = 0. 

E— >oo E— >oo 



This possibility prevents us to claim that for the double well potential (|1.4p with s > s* the two 
branches of asymmetric states concentrated near its two minima, which exist for large E via Theo- 
rem[6j are in fact the continuation of the two branches of asymmetric states emerging from the pitch- 
fork bifurcation along the symmetric branch at E = E* < oo, via Theorem HI Numerical simulations 
strongly support this claim but a proof of Theorem [6] which includes the case x$ = ±oo is needed 
for a rigorous resolution of the problem. Note that, for usual potentials V with \\xa.\ x \_+ o0 V{x) = 0, 
xo = ±oo is actually a degenerate critical point, i.e. 

d k V 

lim —p(x) = Q, k = 1,2,3,.... 

x->±oo ax 

Remark 8. Moreover, in the absence of the spectral hypothesis that L + (ipE, E) has exactly one 
negative eigenvalue, the proof of Theorem part (ii) shows that ue may split in two or more 
functions, depending on the number of negative eigenvalues of L + , which concentrate at points 
further and further apart, i.e there exist u E , u E , . . . , u E 6 H 2 , and y E , y E , ■ ■ ■ , y E 6 ^ such that: 

lim \y> E -y%\=oo for j ^ k, lim \\u E (- - y k E ) - Uoo\\h^ =0, for k = 1,2 ... N, 

E^too E-^oo 

and 

n 

ue=/^u% + ve, where lim \\ve\\li = 0, for all 2 < q < oo. 

k=l 

Such nonlocal bifurcations are excluded for ground states but are relevant for excited states of the 
stationary NLS equation \1.3\) . see fWJjj for partial results in this direction. 
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We now proceed with the proofs of the above theorems. 
Proof of Theorem [3 Recall from (|2.9p : 

£{E)= [ \V^ E (x)\ 2 dx + f V(x)\i; E (x)\ 2 dx + ^— [ \i, E {x)\ 2p+2 dx, 
Jr Jr P + 1 Jr 

and from (pJ0|) 

d£ __ _ E dN 
dE~~ ~dE' 

where N(E) = ||^b||| 2 . Also, from ff27TT|) : 

llV^Hia + / F(x)^ E (x)| 2 c/x + a||^||f+2 2 + S||^||| 2 = 



which can be rewritten: 

8{E) + ^-U E \\% + P l = -EN{E). 
p+1 u 

Differentiating the latter with respect to E we get: 

-^Ue\\%^ 2 = P —N. (4.13) 
dE u — ap 

Now, from ip E a weak solution of (|1.3j) using xip E (x) as a test function and integrating by parts 
we get the Pohozaev type identity: 

-||V^||| 2 + / (V(x)+xV\x))\^ E (x)\ 2 dx + ^—U E \\%+ 2 2 = -EU E \\ 2 L2 (4.14) 

JR P+1 

which added to (12. lip implies: 

2EN+ [ (2V(x) + xV\x))\^ E (x)\ 2 dx = -a^U E \\ 2 L P ^ 2 2 . 

Jr P+1 

The integral term in the identity above can be bounded via Holder inequality: 

- [2||V|| L «» + \\x ■ VV\\ L oo] N< J (2V(x) + xV'(x))\i) E (x)\ 2 dx < [2\\V\\ L °o + \\xV \\ Loo ] N, 

Jr 

hence, for C = 2\\V\\ L : ' + ||xV'||i,°o, we have: 

2(E - C)N(E) < - a ^\\i> E \\f+? 2 < 2(E + C)N(E). (4.15) 

Plugging ()4.15p in (|4.13p and using the notation Q(E) = \ \ip E \\ we get: 

V + 2 Q(E) < dQ < p + 2 Q(E) 



2p E + C ~ dE K ' ~ 2p E-C 
Fix E2 > max{C, E\} and integrate on [E2,E] : 

Q( E i) it? 1 ^1/2+1/p ^ n i-p\ < Q( E 2) lT? ^1/2+1/p 
(E 2 + C)^/p {E + C) - Q(£) -(£ 2 -C)Wp (£ - C) ' 
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Q(E 2 ) < < Q(E 2 ) 



(E 2 + C)V 2 +Vp- EV 2 +Vp ~ (# 2 -<7)V2+1/p' 
which implies (14. 4p since 

QCgg) 

hm = 1. 

(£ 2 +C) 1 /2+l/p 

Now, (I4.5P follows from dividing (14.15P by £' 1 / 2 + 1 / p and passing to the limit E — > 00, while (14.6P 
follows from dividing (|2.1ip by £' 1 / 2 + 1 /p and passing to the limit E — > 00. 

Note that (|4.7p and the fact that t/;^ solves (|1.3p already implies ()4.9|) . Moreover (j4.7|) combined 
with dHU-dUSI) shows that: 

Adding the last two we get (I4.8p . Part (i) is now completely proven. 
For part (ii) we use concentration compactness: 

1. -Vanishing, i.e. lim^^oo ||ue||m = 0, 2 < q < 00, cannot happen. Indeed, assuming the 
contrary, from (|4.9p we get: 



u E = (-A + l)- 1 [-R 2 {V{Rx) - V(0))u E - a\u E \ 2p u E ]. (4.16) 
Hence, using that (—A + : I? 1— > H 2 is unitary, we have: 

\\ue\\h* < R 2 \\V{Rx) - V{Q)\\ L ~\\u E \\v + \<j\\\u E \\ 2 l+l 2 . 
Since the right hand side converges to zero we get ||w.e||#i < ||ub||jj2 E ^f° which contradicts 



2. -Splitting cannot happen. The argument is a slight adaptation of the one we used to exclude 
splitting in the proof of Theorem [3] part (ii) and relies on the hypothesis that L + (^ E , E) has 
exactly one negative eigenvalue for all E G (E\,oo). 

3. -Compactness is the only case possible and implies that for any sequence E n — > 00 there exists 
a subsequence E nk and y^ G R, G H 1 such that: 

lim \\uE k {- - Vk) ~ Uoo\\m = 0. 



As in the compactness part of the proof of Theorem[3]the symmetry of UE k implies that yk G R 
must be a bounded sequence. By possibly choosing a subsequence we have hmfc_> 00 y^ = yoo G 
R, hence 



lim \\u Ek -u OQ (- + yoo)\\m = °- 

k— ¥00 
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By plugging u^ k in (|4.16|) and passing to the limit k — > oo we infer that n QO (- + y^) is an 
even H 1 solution of (|4.ip . hence itoo( - + Uoo) = ^oo given by (|4.2p . 

We have just showed that for each sequence — )• oo there exists a subsequence -E rafe such 
that 

lim - Uoo\\ H \ = 0. 

Since the limit is unique we have: 

lim \\ue - UooWfji = 

and (|4.16|) now shows that the convergence is in H 2 . 

Theorem [5] is now completely proven. □ 
Proof of Theorem^ Since the Frechet derivative with respect to u of G(u, R) : 

D u G(u,R)[v] = L+(u, R)[v] = -v" + R 2 {V{Rx + x ) - V(x ))v + v + (2p + l)a\u\ 2p v 

at u = u^, R = has kernel spanned by u'^ we could use the standard Lyapunov- Schmidt 
decomposition with respect to L + (u oo ,0) to reduce (14, lip to finding the zeroes of a map from R 2 
to M. However the latter will have zero gradient and Hessian at (0,0) because G(0,R) = 0(R 3 ) in 
case (i) and G(0, R) = 0(R 4 ) in case (ii). A generalization of the Morse Lemma and calculations of 
derivatives up to order four will be required to fully analyze the reduced problem. In particular V(x) 
will need to be C 4 or in W 4,co (^i) to be able to carry on the analysis. We avoid this unnecessary 
complication by using a decomposition similar to the one in [7]: 

Lemma 5. There exists e, 5 > such that for any u G L 2 with ||it — Uoo||z,2 < e there exists a 
unique s£R, |s| < S, with the property that: 

u = u 00 (--s) + v, v±u' OQ (--s). 

Moreover the map s is C 2 from L 2 to M. and there exists C > such that 

\s\ < C\\u — Uoo\\ L 2. 



The Lemma follows directly from applying the implicit function theorem to the problem of finding 
the zeroes of the C 2 map F : L 2 x R i— y M given by 

F(u, s) = (u'oq (■ - s), u - «oo (■ - s)) 
in a neighborhood of (uoo; 0), because: 
dF 

— (^,0) = -(u'^, Uoo - Uoo ) + (u'oojU^) = IKoJlia + 0. 

Returning now to the equation (|4.1ip we note that (tioo( - — s), 0), s £ R is the set of all solutions 
when i? = 0. To find other solutions in a if 1 x R neighborhood of (uoo,0) we decompose them 
according to the above lemma: 

u(x) = Uoo(x - s) + v(x), v±u' OQ (--s) (4.17) 
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and rewrite (|4,11|) in the equivalent form: 

G±(v,R,s) = -v" + P ± (s)R 2 [V(Rx + x )-V(xo)}[u 00 {x-s) + v} + v 

+(2p + l)aP ± (s)\ Uoo (x - s)\ 2p v + P ± (s)N(s,v) =0 (4.18) 
G { \(v,R,8) = (u' oc (--s),R 2 [V(Rx + x )-V(x )][u oo (x-s) + v]+N(s,v))=0, (4.19) 

where P±(s) is the projection onto the orthogonal complement of n^ c (- — s) in L 2 , and 

N(s,v) = a\u OD (--s)+v\ 2p (\u O0 (--s)+v)-a\u^(--s)\ 2p u O0 (--s)-(2p+l)a\u O0 (--s)\ 2p v. (4.20) 

Note that, for all x G R and some < i < 1, we have: 

\N(s,v 1 )(x)-N(s,v 2 )(x)\ < (2p + l)2p\a\(\u O0 (x-s)\+m a x{\v 1 (x)\,\v 2 (x)\}) 2p - 1 

x\tvi(x) + (1 - t)v 2 (x)\ \vt(x) - v 2 (x)\, 

hence 

\\N(s, Vl ) - N(s,v 2 )\\ L 2 < (2p+l)2p|cr|(||u 00 || L oo +max{|| Vl || L oo,|| U2 || i oo}) 2p - 1 

x max{||vi||i,oo, ||v 2 [|l,°°}|K - v 2 \\ L 2 

in particular, if we assume ||ui,2||i? 2 < 1 then we can find a constant Cat > such that: 

\\N(s,vi) - N(s,v 2 )\\ L 2 < C , Armax{||t;i|| // 2, ||«2||fl-a}||vi - v 2 \\ H 2 (4.21) 

and 

\\N(s,v 1 )\\ L 2<C N \\v 1 \\ 2 H2 . (4.22) 

Using now the notation L OQ (s)[v] = L + (uoo(- — s),0)[u] = — v" + v + (2p + l)<r|n 00 (- — s)| 2p u we 
can rewrite (14.18P in the fixed point form: 

v = -L^ 1 (s)P ± (s)R 2 [V(Rx + xo)-V(xo)}u 00 (-- S ) 

-L^(s)P ± (s) [R 2 (V(Rx + x )-V(x ))v + N(s,v)] , (4.23) 
= v (s,R) + K s , R (v) 



where: 



L 0O 1 (s) : L 2 n {ttJxj(- — s)}- 1 i y H 2 n {wJxj(- — s)} 1 - is linear, bounded, with bound independent 
on s G R; 

Pj_(s)i? 2 [V%Rx + x )- V^xo)] : # 2 i-> L 2 n{n' 00 (--s)} ± is linear and bounded by 2R 2 \\V\\ L o 
uniformly for s G R; 

P±(s)N(s, •) : H 2 h-> L 2 n{u' 00 (- — s)} -1 is locally Lipschitz with Lipschitz constant independent 
on s G R, Lip(r) < Cjyr in the ball of radius r < 1 centered at origin in i? 2 , see (|4.2ip . 
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Contraction principle can now be applied to (|4.23|) in the ball: B(0,r) = {v £ H 2 : \\v\\jj2 < r} 
provided ||^o(s, R)\\h 2 — r /2 and K s ^r is a contraction on -6(0, r) with Lipschitz constant Lip < 1/2. 
Based on the above estimates it suffices to require r < 1 and: 



2|. ]|. 

2-R ||A>o IU 2 i-KH" 2 l|>^ 1 1 -L°° I Poo 1 1 if 2 < - 

1 



\\L^\\ L 2^ H 2(2R 2 \\V\\ L ~+C N r) < | 
which can be accomplished by choosing: 

r = min{l, (2\\L^ ll^^j^Cjv+l/lltioc lli^) -1 }, R < V^ll^^lli^^ll^llL^ll^oollj^) -1 = -Ro 

Hence (|4.23p . and consequently (|4.18[) . has a unique solution v = v(s,R) in B(0,r) for each sgK 
and < R < Rq. This solution depends C 1 on s and R since (|4,23[) is C 1 in these parameters. 
v(s,R) can be obtained by successive approximations: 

H 2 

v (s,R), vi(s,R) = v (s,R) + K SiR (v (s,R)), v 2 (s,R) = v (s,R) + K SjR (vi(s,R)), . . . -)■ u(s,i?) 
Moreover, from the contraction principle we have: 

\\v(s,R) - vi(s,R)\\ H 2 < LiP \\vi - vq\\ H 2 < ||L^ 1 || L 2 l _ >H2 (2 J R 2 ||y|| L oo + C n \\v \\h 2 )\\vo\\h 2 - 

1 — Lip 

In what follows we will show that for \s\ < 1 we have: 

M , mil f <3(i? 3 ) if V'(x ) + 

lh(,,i?)||H2 = { \ Ri > jf^jlo (4-24) 
where in the second case we assume y is twice differentiable at xq. Consequently: 

=V^) + { ^ ify'U) = o ^ 
For (|4.24|) we use the definition (|4.23|) of v : 

\\voWm < R 2 \\L^\\ L ^m\\[V(Rx + x ) - V(x )]u oo (- - s)\\ L 2, 

and 

\\[V(Rx + xo)-V(x )}u 00 (--s)\\ 2 L 2 = [[V(Rx + x )-V(x )] 2 u 2 oo (x-s)dx 



/ u 2 00 (x){V(Rx + Rs + xq) - V(x )] 2 dx. 
Jr 



Since V{x) is differentiable at x = xq and bounded for almost all x, there exists C\ > such that 
almost everywhere: | V(a;) — V(xo)| < C\\x — xq\. Moreover if V(x) is twice differentiable at x = xq 
with V'(xo) = then there exists C 2 > such that almost everywhere: \V(x) — V(xq)\ < C 2 \x— xo| 2 . 
Plugging in above and using J" R Ix^-u^x)^ < 00, n = 0, 1, 2. . . due to the well known exponential 
decay of Uqo, we get (|4.24[) . 
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Since (|4.1ip is equivalent with the system (|4.18|) - (|4.19p . and (|4.18j) has the unique solution v(s, R) 
in a neighborhood of zero for all s € R, < R < Rq, all solutions of (|4.1ip in a small neighborhood 
of (uoo,R = 0) will be given by the solutions of (|4.19|) with v = v(s,R) and (s,R) in a small 
neighborhood of (0, 0) : 

G l{ (s,R) = G ll (v(s,R),R,s) = (u' 00 (--s),R 2 [V(Rx + xo)-V(x )}u 00 (--s)) (4.26) 
+ « (- - s),R 2 [V(Rx + x ) - V(x )]v(s, R) + N{s, v(s, R))) = 0. 

Note that by the estimate f04~|) - (j4T2"5|) we have, for V'(x ) / and \s\ < 1 : 

(n' 00 (- - s),R 2 [V(Rx + so) - V(x )]v(s, R) + N(s, v(s, R))) = 0(R 6 ) 

while 

(u' oo (--s),R 2 [V(Rx + x )-V(x )]u oo (--s)) = R 2 [ vf x (x)u 0O (x)\V(Rx + Rs + x ) -V(xo)]dx 



/ u 2 X) {x)V' {Rx + Rs + xo)dx 
Jr 



R*_ 
2 

n^)hoo\\i, R3+o{R3) 



where in the last step we used lim^o V'(Rx + Rs + xq) = V'(xq) and Lebesque dominated con- 
vergence theorem. Hence, for V'{xq) ^ and |s| < 1 : 



^ {l (s,R) = - V ' {X0) ^R^oiR 



°°llL 2 r>3 



which has only the R = 0, s £ R in a small neighborhood of (s, R) = (0, 0). However, for V'(xq) = 
and |s| < 1, we have 

R? f 

(u' 00 (--s),R 2 [V(Rx + xo)-V(xo)}u 00 (--s)} = -— / u^^V (Rx + Rs + x )dx 

1 Jr 

R A [ u 2 OD (x)(x + s)dx + o(R 4 ) 
Jr 



V"(x ) 



d4 V"(x ) 2 4 

= ~R \\Uoo\\ L 2S + o{R ), 

hence 



G {l (s,R) = R- A G\ { (s,R) = -^M|| Uoo ||2 2S + (i) = o 



for which the implicit function theorem can be applied at (s,R) = (0,0) where d s Gii(0, 0) = 

— V 2^°^ ||t/.oo ||^a 0- Note that the o(l) term is differentiable with respect to s because v(s,R) is 
and it remains o(l) in R. 

The estimate 

IMlls = llwoollia + " i) V"(x )R 4 \\x Uoo (x)\\ 2 L2 +o(R 4 ). 
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can be obtained from ur = n 0O (- — s(R)) + v(s(R), R), see (|4.17|) . which implies: 

\\ u r\\l 2 = ll n oo|li2 + 2(uoo(- - s),v) + \\v\\ 2 L 2 = H^oo Hi.2 + 2(1400 (• - s),v ) + 0(R 6 ) 
From (|4.23p we have: 

(«oo(- - s),v ) = -(L c ^ 1 ti oc (- - s),R 2 [V(Rx + xq) - V(xq)]u 00 (- - s)) 
= -W V "( x o) [ ( x + s) 2 Uoo (x)L^[ Uoo ]dx + o(i? 4 ), 

1 JR 

and the integral term can be computed using — L^Uqo] = 2^«oo + ^xu'^. 

The main part of the theorem is now finished. It remains to prove the spectral properties of the 
linear operator: 

L+(R)[<f>] = D u G(u(R),R)[(j>} = -f + R 2 (V(Rx + x ) - V(x ))<f> + + (2p + l)<r|u(i2)| 2p <£, 

where u(R) = Uoo(- — s(R)) + v(s, R). It is well known that: 

L+(Q)[<f>] = Loo(O)[0] = -0" + </» + (2p + l)a| Uoo | 2p 0, 

has exactly one, strictly negative eigenvalue which is simple and zero is the next eigenvalue which 
is also simple with corresponding eigenvector u'^. Since u(R) depends continuously on R, L + (R) is 
continuous with respect to R in the resolvent sense, hence the isolated eigenvalues and eigenvectors 
depend continuously on R. To establish the sign of the second eigenvalue we use the following 
expansion for its eigenvector: 

<f>{R) = (1 + aiR^uU- ~ s(R)) + $(R), $(R) 1 n' 00 (- - s(R)), 
and the eigenvalue equation: 

L+(R)[<f>(R)] = \(R)cj ) (R) 

which is equivalent to: 

X(R)4>(R) = L 00 (s)4>(R) 

+ [R 2 (V{Rx + x ) - V(x )) + (2p + l)2pa\ Uoo {- - s)| 2p -y S , R)] 4>{R) + DN(R)4*{R) 

V * ' 

y(s,i?)=o(R 4 ) 

where 

DN(R) = (2p + l)a {\ Uoo (- -s) + v(s,R)\ 2p - \ Uoo (- - s)| 2p - 2p\ Uoo (- - stf^v&R)) 

( v(s,R) if p = 1/2 
= < 0(R 8p ) if 1/2 < < 1 
{ 0{R S ) iip> 1. 

We project the eigenvalue equation onto u^O — s (^)) an d its orthogonal complement in L 2 to 
obtain the equivalent system of two equations: 

- \{R) + P±(s)F(s, fl) + P ± (s)DN(R)) $(R) 
= ~P±(s) [V(s, R) + DN(R)] (1 + a(R)) Uoo (- - s) 
A(i?)(l + a(i?))|KJ|| 2 = (u> 00 (--s(R)),V(s,R) ( f>(R))+o(R 4 ) 
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As before, from the first equation we deduce <&(R) = 0(R 4 ), while replacing now in the second one 
<j>{R) = (1 + o(l)K 00 (- - s{R)) + 0(i? 4 ), we get: 

\{R) = IKoH^Voot- - S ),^^Koo(- - *)> + o(i? 4 ). 

To calculate 

(^(•-^(s^K^-s)) = (n / oo (-- S ),i? 2 [F( J Rx + x )-y(x )]n' oo (.- S )) (4.27) 

+ (n' 00 (- - a), (2p + l)2H«oc(- - s^V^-RKoG - «)) 

we use i>(s, i?) = t>o(s, R) + 0(R 6 ) and the equation satisfied by t>o, see f|4.23[) : 

Loo(s)v = -P ± (s)R 2 [V(Rx + x ) - V(x )] Uoo (- - s). 

Taking its space derivative: 

L^sK + (2p + l)2H«oo(- " s)| 2p-1 <o(- - «H 
= -R 3 V'(Rx + x )noo(- - s) - fl 2 [V(ita + x ) - F(x )]u' oo (- - s) 
+(u' 00 (- - s),R 2 [V(Rx + x ) - V(x )] Uoo (- - s)X(- - s) 

then its scalar product with w^ 50 (- — s) we get 

(<4(- - a), (2p + l)2jw|«oo(- - s)| 2p -%)( S , - s)) 

= -i? 3 «o0 - s),V'{Rx + x ) Uoo (- - s)) - i? 2 (n , 00 (- - s) 5 F(ita + xq) - ^(x )]n' oo (- - s)) 

which plugged into (|4.27|) leads to: 



(uU- ~ s),V(s, R)uU- ~ a)) = -R 3 (uU- ~ s), V'(Rx + x )n oo (- - a)) + 0(R 6 ) 

2 



^V"(x )\\ Uoo \\l 2 +o(R 4 ) 



All in all we have: 

\(R) = ~V"(x )R 4 + o(R% 

which shows that the second eigenvalue of L + {R) becomes negative (respectively positive) if x$ is 
a local maxima (respectively local minima) for the potential V(x). The theorem is now completely 
proven. □ 



5 Numerical results 



To illustrate the results and further investigate the behavior of the ground state branches we have 
performed a series of numerical computations on the equation 

+ -V(x)i(>(x) - \ip(x)\ 2p 4>(x) + Ei(>(x) = 0, 

which is equivalent with f 1 1 . 3 [) with o = — 2 if one regards the eigenvalue parameter E in this section 
as being half the parameter used in the previous sections. The potential V = V s is given by (jl.4p . 
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When s < ~ 0.6585, the potential is a single well but it becomes a double well for s > s*. We also 
recall that p = p* ~ 3.3028 is the theoretical threshold power for the nonlinearity that separates 
different pitchfork bifurcations in the regime s — > oo, see Corollary [21 The principal conclusions of 
our investigations for p = 1, p = 3, and p = 5 can be summarized as follows: 

1. When the potential is a single well (that is s < s*) ; the symmetric ground state exists for all 
E > Eq and the operator L + at this states has a single negative eigenvalue for all E > Eq. 
There are no bifurcations along this branch which is consistent with the result in 

2. When the potential is a double well (that is s > s*), the symmetric ground state exists for 
all E > Eq and the second negative eigenvalue of L+ along the branch of symmetric states 
emerges for E > E*, where E* depends on s and E* > Eq. The asymmetric states bifurcate 
at E = E* and exist for all E > E*. The second eigenvalue of L+ along the branch of 
asymmetric states is positive for all E > E*. The numerical results are in agreement with 
Theorem [TJ Furthermore they show that there are no other bifurcations along these branches 
past £7*, that, as E — > oo, one branch of asymmetric states localizes in the left well while 
the other localizes in the right well, and, modulo re-scaling, they both converge to the NLS 
soliton localized in the left, respectively right, minima of the potential, see Remark [6l 

3. When p < p*, the pitchfork bifurcation is supercritical and the branch of asymmetric states 
has bigger I? norm than the one for the symmetric state at E = E*. When p > p*, the 
pitchfork bifurcation is subcritical and the branch of asymmetric states has smaller L 2 norm 
for E ~Sl E* than the one for the symmetric state at E = E*. The numerical results are 
consistent with Corollary [2] but also suggest that the separation between wells does not have 
to be large, i.e., as soon as the potential has two wells, the supercritical/subcritical character 
of the bifurcation is controlled by the nonlinearity. 

The conclusions are showcased in the following five figures. Figure Q] illustrates the cases of 
s = 0.6 and s = 0.7, for p = 1, that straddle the critical point s* « 0.6585. The top panel presents 
the dependence of the squared 1? norm of the symmetric state on E, while the bottom panel 
presents the second eigenvalue of the operator L+ as a function of E. It is clear that for s < s*, the 
second eigenvalue of L + tends asymptotically to 0, without ever crossing over to negative values 
(solid line), while for s > s#, such a crossing exists (dashed line), occurring for E* 10.68. On 
the other hand, to examine whether a secondary crossing may exist for larger values of E, we 
have continued the s = 0.7 branch to considerably higher values of E in the bottom right panel of 
the figure, observing the eventual convergence of the eigenvalue to A = 0, without any trace of a 
secondary crossing to positive values. 

Figure [2] further clarifies the bifurcation structure of the asymmetric states for s = 0.7>s*, p = 
1. As a relevant diagnostic, we monitor the location of the center of mass of the solution 

f R x\iiE\ 2 dx 

XCM = rr~i \ia ■ 

J R WE\ 2 dx 

We can clearly see from the top left panel that beyond the critical threshold of E* ~ 10.68, two 
asymmetric states with xqm ^ 0, corresponding to a < and a > in Theorem [H bifurcate out 
of the symmetric state with xcm = 0, with the latter becoming unstable as per the crossing of 
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Figure 1: The top panel shows the dependence of the squared L 2 norm of the symmetric state on 
the parameter E for p = 1. The bottom left panel shows the trajectory of the second eigenvalue of 
L + for the cases of s = 0.6 < s* (blue solid line) and s = 0.7 > s* (red dashed line). The bottom 
right panel shows an expanded plot of the second case for considerably larger values of E. 
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Figure 2: The top left panel shows the pitchfork bifurcation of asymmetric states for s = 0.7 > s* 
and p = 1. Two asymmetric (xcm ^ 0) states emerge for E > E* « 10.68. The second eigenvalue 
of L_|_ for such asymmetric states is positive as shown in the top right panel. The bottom panel 
shows symmetric and asymmetric states versus x for E = 15. 

the second eigenvalue of L + to the negative values. For the asymmetric branches with xcm ^ 
emerging past the bifurcation point, the second eigenvalue of L + is shown in the top right panel of 
the figure, with its positivity indicating the stability of asymmetric states. These panels corroborate 
the supercritical pitchfork bifurcation scenario Q > 0, R > in Theorem HJ The bottom panel 
shows both symmetric and asymmetric states for E = 15. 

The relevant computations are repeated for higher values of p. The corresponding numerical 
results for p = 3 are shown in Figure O Again illustrating the cases s = 0.6 < s* and s = 0.7 > s*, 
we observe that a crossing of the relevant eigenvalue occurs in the latter but not in the former. 
Notice that in the latter case of s = 0.7, as shown in the bottom left panel of Fig. O the second 
eigenvalue crossing occurs for E* ~ 7, i.e., for a smaller value of E than in the p = 1 case. 
Generally, we have found that the higher the p, the earlier the relevant crossing occurs and also the 
more computationally demanding the relevant numerical problem becomes, as the solution narrows 
and it becomes challenging to appropriately resolve it even with a fairly fine spatial grid for large 
values of E. This is clearly illustrated in the bottom right panel of Fig. [3l where it can be seen that 
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Figure 3: The top and bottom left panels are similar to Fig. [IJ but now for p = 3. The blue solid 
line corresponds to the case of s = 0.6 < s*, while the red dashed to s = 0.7 > s*. The bottom right 
panel shows the importance of sufficiently fine discretization in resolving this second eigenvalue for 
large E. Here the red dashed line corresponds to a spatial grid spacing of Ax = 0.025, while the 
solid blue line is obtained for Ax = 0.0125. 

in the absence of sufficient resolution (dashed line for larger spacing of the spatial discretization), 
a spurious secondary crossing is observed for the second eigenvalue of L + . This secondary crossing 
is eliminated by finer discretizations (solid line). 

It should be pointed out that the symmetric states become unstable in the supercritical case p = 3 
for E w 0.95, s = 0.6 and E r< 0.85, s = 0.7, due to the change of the slope of N(E) = \\i/je\\ 2 L 2 (see 
the top panel of Fig. [3]). Therefore, the asymmetric states bifurcating from the symmetric ones at 
E* w 7 will also be orbitally unstable, because Q > 0, \'(E*) < 0, and N'(E*) < will imply R < 
in Theorem 01 However, as s becomes larger, the value of E* becomes smaller and E* — > Eq as 
s — > oo, see Corollary [2j Figure H] shows the dependence of the squared 1? norm for the symmetric, 
asymmetric, and anti-symmetric stationary states for p = 3 and s = 10. The pitchfork bifurcation 
occurs while the slope of HV'eII^ is still positive and leading to orbitally stable asymmetric states 
and a change of stability along the symmetric states from stable for E < E* to unstable for E > E* 
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Figure 4: The graph shows the dependence of the squared I? norm of the stationary states on the 
parameter E for p = 3 and s = 10. The left most (blue solid) line are the anti-symmetric (exited) 
states emerging from zero at the second lowest eigenvalue E\ = 0.191042. Almost on top of it is 
the (blue solid) line of symmetric ground states emerging from zero at Eq = 0.191046. The latter 
bifurcates at E* <J Eq into the dashed red line (symmetric states) and the right most blue solid line 
(asymmetric states). At E ~ 0.5908 the slope of the L 2 norm becomes negative for all branches. The 
solid lines denote linearly stable branches, while dashed ones denote unstable branches. 

as stipulated in Corollary [2j Note that the numerical simulations suggest that the asymmetric 
branches can be continued for all E > E*, hence the slope of their L 2 norm square will change for 
large E and they will become unstable, see Fig. [Hand Remark (6[ 

The numerical results for p = 5 > p* and s = 4 > are shown in Figure [5j In this case, the 
subcritical pitchfork bifurcation occurs at E* « 0.196 at the positive slope of ||V , _b||^2 with respect 
to E. The top panels show the behavior of squared 1? norms for symmetric, anti-symmetric, and 
asymmetric stationary states. The blowup on the top right panel shows that the asymmetric states 
have decreasing 1? norm for E* < E < E* , corresponding to the case Q > 0, R < in Theorem 
IH and in agreement with Corollary [21 However, their L 2 norm becomes increasing for E > E*, 
where E* ~ 0.202, hence the asymmetric states are orbitally stable in this regime. The bottom 
panels show the squared eigenvalues of the stability problem associated with the symmetric (left) 
and asymmetric (right) states. The symmetric state is unstable for any E > E* (because of the 
second negative eigenvalue of the operator L + ). It becomes even more unstable for E > E*, where 
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Figure 5: The top panels show dependence of the squared L 2 norm of the stationary states on the 
parameter E for p = 5, s = 4, where the right panel is a blowup of the left panel. The leftmost 
branch is the anti-symmetric (excited) states and the other branch corresponds to the symmetric 
ground states which bifurcate into asymmetric states. The bottom panels show the squared eigen- 
value of the linearization spectrum associated with the symmetric (left) and asymmetric (right) 
branches. The insert on the bottom right panel gives a blowup of the figure to illustrate instability 
of asymmetric states near the subcritical pitchfork bifurcation. The solid lines on the top panels 
denote linearly stable branches, while dashed ones denote unstable branches. 

E* tss 0.32, when another unstable eigenvalue appears (because of the negative slope of the 1? 
norm). The asymmetric state is unstable for E* < E < E* (because of the negative slope of the I? 
norm) and stable for E > E* . We note that the asymmetric state becomes unstable past E ~ 0.32 
because of the negative slope of the L 2 norm, similarly to the branch of symmetric states and 
consistent with Remark [6l 

While the results presented herein formulate a relatively comprehensive picture of the one- 
dimensional phenomenology in the context of double wells, some questions still remain open for 
future investigations. One of them, raised in Remark [TJ is associated with states of the form of Uoo 
spatially concentrated at ±oo. Our numerics for the potential (jl.4p has not revealed such states 
presently, but it would be relevant to lend this subject separate consideration. Additionally, we have 
not seen the case Q < in Theorem 01 i.e. the eigenvalues E are decreasing along the asymmetric 
branch. Another important question concerns the generalization of the results presented herein 
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to higher dimensional settings. There, the bifurcation picture is expected to be more complicated 
when the eigenvalues crossing zero are not simple. 
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